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Abstract 

An analysis of current computational fluid dynamics capabilities in predicting mean lift forces 
for two dimensional foils is conducted. It is shown that both integral boundary layer theory 
and Reynolds Averaged Navier Stokes algorithms provide the same over-prediction of lift forces 
when properly converged. It is also shown that the over-prediction is insensitive to turbulence 
model details. 

Experimentation and computational fluid dynamics modeling show that discrete vortices 
are shed with significant sizes and distinct frequencies. These vortices are shown to result in 
significant cfd prediction errors when they are asymmetric in size or shape. Inaccuracies in flow 
predictions in the near wake appear to result in an effective change in the Kutta Condition due to 
pressure biasing associated with vortex asymmetry. The net result is a consistent overprediction 
of mean lift. 

Based on an analysis of over 1000 historical experiments an empirical model is developed to 
allow the error in predicted lift coefficient to be anticipated based on the local flow conditions 
at the trailing edge of the foil. A series of experiments are conducted and reported to test the 
accuracy of the empirical model. The result is a significant improvement in mean lift prediction 
and pressure profile for both RANS and IBLT. 

Thesis Supervisor: David Burke 

Title: Senior Lecturer, Department of Ocean Engineering 
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Chapter 1 


Introduction 


The thrust produced by a marine propeller is a result of a propeller blade “wing” section 
oriented such that the relative motion caused by rotation produces “lift” forces in the direction 
of ship motion. It is for this reason that much attention has been given by the marine industry 
to analysis and prediction of two dimensional foil section performance. 


1.1 History 

1.1.1 Efforts Prior to Computer Modeling 

A great deal of effort has been given to wing section development by the aircraft industry. The 
early stages are described in Reference [1]. The quote here is from 1949. 

“Until recently the development of wing sections has been almost entirely empirical. 

Very early tests indicated the desirability of a rounded leading edge and of a sharp 
trailing edge. The demand for improved wings for early airplanes and the lack of 
any generally accepted wing theory led to tests of large numbers of wings with shapes 
gradually improving as the result of experience. The Eiffel and early RAF series 
were outstanding examples of this approach to the problem. 

The gradual development of wing theory tended to isolate the wing-section problem 
from the effects of plan form and led to a more systematic experimental approach. 

The tests made at Gottingen during the First World War contributed much to the 
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development of modern types of wing sections. Up to about the Second World War, 
most wing sections in common use were derived from more or less direct extensions 
of the work at Gottingen . During this period, many families of wing sections were 
tested in the laboratories of various countries, but the work of the NACA was out¬ 
standing. The NACA investigations were further systematized by the separation 
of the effects of camber and thickness distribution, and the experimental work was 
performed at higher Reynolds numbers than were generally obtained elsewhere. The 
wing sections now in common use are either NACA sections or have been strongly 
influenced by the NACA investigations.” 

The wing theory discussed in this quote is strictly inviscid. The results of these investiga¬ 
tions are very well documented and freely available. In the process of conducting this research I 
have considered many experimental data points and the characteristics of several facilities. Af¬ 
ter having reviewed the experimental procedures, opinions of others, and the facilities utilized, 
it is my opinion that the most accurate experimental two dimensional wing section charac¬ 
teristics available are those conducted in the NACA two-dimensional low-turbulence pressure 
tunnel. The majority of the wings tested in this facility are presented in Reference [1] and a 
few more are presented in NACA Technical Notes and Technical Reports published after this 
reference. 

1.1.2 Exploring Capabilities in Computer Modeling 

As computer speed and software progress the viscous effects are increasingly being considered 
in a general way. The time averaged effects of viscosity can be accounted for by modifying 
the local boundary conditions at the edges of the inviscid domain according to the integrated 
effect of the boundary layer or by turbulence modeling throughout the domain. A couple 
examples where the former method has been successful are PBD [17] and XFOIL [8] . The 
latter method is typically associated with RANS solvers such as DTNS [13], UNCLE [3] and 
CFDSHIP [23]. 

Because of the circumferential non-uniformity a marine propeller typically operates in, a 
study was conducted at MIT concerning the need for modifications to the Kutta Condition 
and other flow models in the presence of unsteady inflow conditions [22]. Experiments were 
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conducted and blind predictions performed by the authors of a number of flow codes. The panel 
of experts who modeled this experiment in the computational flow domain found a mixture of 
encouraging and discouraging results. They found that in general, it was possible to model the 
unsteady behavior and achieve results consistent with the variations observed in the experiment. 
On the other hand, they found that the zeroth harmonic (or mean behavior) showed significant 
deviation when typical turbulence models were used. One of the publications following this 
project put it this way: [23] 

The zeroth harmonic, which was the same as steady flow calculations, displayed typ¬ 
ical foil response and showed that turbulence model corrections for pressure gradient 
were necessary for accurate solution. 

and 

...This level of agreement, which required pressure-gradient modifications to the tur¬ 
bulence model, appears to be a general assessment of the current capabilities of 
isotropic turbulence models since it is consistent with the other participants m the 
FFX workshop and overall results found in the literature. 

1.1.3 Efforts to Reconcile Experiments and Computations 

As a result of these inconsistencies a study of the steady behavior of two dimensional foils and 
the associated computations was undertaken at the MIT water tunnel. A foil with a more 
drastic cup than typical for the type of anti-singing trailing edge in the flapping foil experiment 
was designed and experimentally evaluated. Blind tests were conducted by the authors of the 
same computational fluid dynamics codes. In an engineer’s thesis by John Dannecker [7] the 
results of the different flow codes were shown to have major variations between flow codes. 
Also, each of the codes over-predicted the experimental results by between 25% and 100% for 
an experimental lift coefficient of about 0.5. 

The radical errors and consistent over prediction by the flow codes motivated a more exten¬ 
sive validation of the experiment. The questions about the experiment and associated steady 
predictions were the sole topic of a masters thesis by Jergen Jorde [15]. He concluded that the 
original experiment was correct and the most probable cause for the discrepancy was incorrect 
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Figure 1-1: Vortex Shedding Patterns Behind Several Foils at High Reynolds Number 


modeling of tunnel walls, but there was no final closure on the differences between the codes 
and the over prediction of lift forces. 

In parallel with the work of Jorde, Kimball and Shearer conducted a study of foil trailing 
edge details as it pertains to cavitation performance improvement [19]. Their study was 
focused on cavitation performance design and experimental validation, so the lift predictions 
were not explicitly compared to the experiment. A part of this experiment included high speed 
photography which revealed the wake structure shown in Figure 1-1. It was interesting to note 
that for the blunt and beveled trailing edge foils there was a distinct repeated organization to 
the vortex shedding in the wake. It appears that there is significant viscous-inviscid interaction 
at a frequency that would make accurate predictions with any sort of time averaged flow code 
questionable. 

Following these efforts Fairman [11] conducted a survey of the experimental and computa¬ 
tional efforts on various foil sections evaluated in the MIT water tunnel. He updated both 
experimental and computational efforts for the foils discussed in the previous history at MIT. It 
was shown that the three dimensional effects at the tunnel walls do not effect the ability of the 
foil to represent a two-dimensional space near the mid-span. He also updated the computational 
results to fully include the effects of the upper and lower tunnel walls. 







1.2 Organization 

1.2.1 Overview 

The current research attempts to bring to closure many of the issues surrounding the work on 
two dimensional hydrofoils at the MIT Water Tunnel. The following summary of each chapter 
is provided to allow the reader to navigate to the section of personal interest: 

1.2.2 Chapter 2 

A brief review of the physics associated with the various forms of computational fluid dynamics 
in use today is conducted. The governing equations and methods of implementation are briefly 
discussed with references given for further evaluation. 

1.2.3 Chapter 3 

A comparison between experiments performed by other facilities and computational predictions 
by the integral boundary layer theory code MSES is performed. The quality of the historical 
experiments is evaluated and appropriate references given to allow further evaluation. A 
parametric evaluation is conducted and a simple parametric model is proposed to account for 
differences between experimental and computational lift. 

1.2.4 Chapter 4 

An experimental validation of the parametric model suggested in Chapter 3 is performed. It 
is shown that the full viscous solution combined with the parametric fit provides consistently 
outstanding results. 

1.2.5 Chapter 5 

It is shown that Integral Boundary Layer Theory and Reynolds Averaged Navier Stokes (RANS) 
codes provide the same over prediction in lift coefficient when fully converged. It is shown that 
several RANS turbulence models also provide the same over prediction of lift coefficient. 
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1.2.6 Chapter 6 


Possible assumptions common to these codes that could potentially be in error are considered. 
It is shown that the most likely erroneous assumptions can be related to known coherent vortex 
shedding at the foil trailing edge. 

1.2.7 Chapter 7 

An analysis of the effects of vortex shedding and unsteadiness on time average lift coefficient 
is provided. It is shown that it is possible to provide a complete time accurate calculation of 
the flow domain using unsteady RANS. The average of the time accurate calculation results 
in lift predictions which are consistent with time average measurements. The conclusion is 
drawn that the need for a parametric adjustment to the full viscous solution is caused by the 
fact that the vortex shedding patterns represent an unmodeled physical phenomenon relevant 
to determining time averaged lift. 

1.2.8 Chapter 8 

A brief summary of the conclusions and how they can be practically utilized is presented. 
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Chapter 2 


Overview of Turbulence Modeling 
for Computational Fluid Dynamics 


Vast quantities of books and papers are available on the topic of turbulence modeling. When 
this research is presented it is common that the response of the turbulence modeling expert is 
to wonder if the new turbulence model that they are working on will solve the entire problem. 
Certainly this may someday be true, but for the time being this research is using turbulence 
models which are widely available for engineering applications. Clearly much progress has been 
made in turbulence modeling and many types of flows seem to have been mastered. 

The Reynolds number at which marine propellers typically operate is very large; therefore, 
there is clearly a need for turbulence modeling. The manner in which turbulence models 
determine the Reynolds-stresses varies, but all of them essentially use a model for turbulence 
behavior coupled with experimentally determined constants. Reference [14] is a total of 34 
pages and is an excellent summary of turbulence models in current use. 


2.1 Integral Boundary Layer Models 

The integral momentum equation (Equation 2.1) was first derived by Karman (1921) and is 
taken here from White. [29] 
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( 2 . 1 ) 



This is the most simple and robust type of turbulence model. This method does not 
explicitly consider the changing conditions across the boundary, but makes use of the shape 
factors and integral quantities to represent turbulence quantities. This allows the partial 
differential equations needed to explicitly solve the entire flow domain to be replaced by a set 
of ordinary differential equations representing the integrated quantities. 

This type of turbulence model can then be coupled with an Euler Solver (as in MSES [10] 
) or a Vortex Lattice Method (as in PBD [17] and XFOIL [9]) to provide very fast and stable 
solutions. It is shown in this research, that in the case of two-dimensional foils, these models 
are every bit as accurate as more general eddy viscosity models and require a small fraction of 
the time and effort to compute a flow solution. (Chapter ??) 


2.2 Eddy-Viscosity Models (EVM) 

This type of turbulence model assumes a single scalar quantity (eddy viscosity) coupled with 
velocity gradients can be used to adequately represent the Reynolds-stresses throughout a flow 
field (Equation 2.2). Within this concept is the notion that turbulent transport and molecular 
transport contain an analogy useful for solving the flow properties, but eddy viscosity is NOT 
a fluid property and varies a great deal throughout the flow domain. Most modern Reynolds 
Averaged Navier Stokes (RANS) codes use this type of turbulence model. 
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( 2 . 2 ) 


-v!-v f - 

i 3 




Vt = 


dui duj 
Ut ^ dxj dxi 
Reynolds stress 

eddy viscosity 



k = turbulent kinetic energy 


The manner in which turbulence models in this category determine the eddy viscosity to be 
used is varied. In essence, the dimensions of eddy viscosity and understanding of the turbulence 
problem lead to the following non-dimensionalization found in many references such as [20]: 


Vt 

f(x,y ) 
V 
L 


f(x, y)xV x L 

dimensionless spatially varying function 
turbulent velocity scale 
turbulent length scale 


(2.3) 


The following is a summary of how many modern turbulence models estimate the parameters 
involved in the eddy viscosity turbulence models. 

• algebraic models 

These models use purely algebraic relationships to define all three of the constants in 
equation 2.3. Examples of these types of models are: 

— mixing length model 
— Cebeci &; Smith model 
— Baldwin & Lomax model 


• one-equation models 

These models use one differential equation to model some characteristic property of the 
turbulence (typically eddy viscosity or turbulent kinetic energy) and then supplement 
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with algebraic equations for the remainder of the constants. Examples of this type of 
model include: 

- Spalart and Allmaris 

- k-1 model of Bradshaw et al 

• two-equation models 

These models use two full sets of differential equations including some sort of production, 
loss and diffusion of eddy viscosity to model the eddy viscosity throughout the domain. 
Examples of this type of model include: 

- k-e Model (Differential equations are solved for turbulent kinetic energy and dissi- 
pat ion.) 

k-u; Model (Differential equations are solved for turbulent kinetic energy and turbu¬ 
lence dissipation rate.) 

2.3 Other approaches 

Differential Reynolds Stress Models (DRSM) solve a complete set of differential equations for 
each component of the Reynolds Stress Tensor. In two dimensions this approach requires three 
sets of differential equations. In three dimensions six sets of differential equations are required. 

Large-Eddy Simulations (LES) model only small scale turbulent structures and large scale 
structures are calculated based on an unsteady time integration. 

Direct Numerical Simulations solve the unsteady Navier-Stokes equations without applying 

a turbulence model. Very small time steps and extremely fine grids are required for this 
approach. 

2.4 Turbulence Model Limitations 

The turbulence models that are readily available for marine propulsion application utilize either 
the integral boundary layer or eddy viscosity approach. The following limitations exist for 


18 


Flow feature that can lead to in¬ 
accurate prediction[14] 

Potential Implication for Pre¬ 
dicting Two Dimensional Lift 

flows with abrupt changes in the ve¬ 
locity field 

trailing edge pressure side where ex¬ 
tremely thin boundary layer exists 

flow along strongly curved surfaces 

trailing edge of a foil with separation 

flow along strongly curved surfaces 

trailing edge of a foil with blunt trail¬ 
ing edge 

flow along strongly curved surfaces 

leading edge of a foil 

stagnant flow regions 

leading edge and trailing edge 

low frequency unsteadiness 

trailing edge in the presence of vortex 
shedding 


Table 2.1: Turbulence Model Weaknesses 


existing one equation and two equation eddy viscosity models and integral boundary layer 
turbulence models: 

• The turbulent viscosity is assumed to be a scalar quantity. 

• Turbulence anisotropy is not explicitly taken into account. 

• A linear relation exists between the Reynolds stresses and the gradient of the mean velocity 
field. 

Many sources, such as [14] and [4] discuss some of the particular weaknesses of these types 
of turbulence modeling. In particular, these restrictions can lead to inaccurate prediction of 
flows in regions which exhibit specific types of flow characteristics. In the first column of Table 
2.1 is a list of flow features Henkes [14] specifically states result in the potential for inaccurate 
predictions. In the second column of Table 2.1 is an analysis of how this type of flow might 
effect the prediction of lift on a two-dimensional foil. 

These weaknesses are combatted using wall functions and other creative turbulence model 
modifications. Often these modifications are specific to an application or a location in the 
domain. Many publications and books have sought to provide closure to these issues in a 
myriad of ways. An expert in turbulence modeling, Henkes summarizes his view of the status 
of turbulence modeling in a book published as recently as 1998 [14]: 
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As the Reynolds number for the flow around the aerofoils and fuselage of aircraft is 
usually very large, most of the boundary-layer flow along the surface will be turbulent. 
It is not possible to compute all the details of the turbulent motion within a reasonable 
turn-around time (say several hours); therefore for all practical computations the use 
of turbulence models is mandatory. 

Up to now most designers and developers of aerofoils compute turbulent flows with 
older generation models, like the algebraic models of Cebeci-Smith, Baldwin-Lomax 
or Johnson-King. More recently some industries and research institutes (like NLR) 
have also applied two-equation models. All these models assume the existence of 
a single turbulent velocity and length scale, which is only approximately right for 
so-called equilibrium flows, such as attached boundary layers in a zero or moderate 
streamwise pressure gradient. Indeed those models have shown to be rather inaccu¬ 
rate for non-equilibrium flows, like separating boundary layers occurring in high-lift 
configurations; examples of such configurations are aerofoils under high angle of at¬ 
tack, and multi-element aerofoils. Therefore all these models are only accurate in a 
limited number of flow types, and they do not meet the high-accuracy requirements 
for a unde range of configurations. As a result of this, experiments in wind tunnels 
still serve as the major source of design information. 
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Chapter 3 


Integral Boundary Layer Theory 
Predictions Compared to NACA 
6-Series Experimental Data 


The experiments referenced here are fully reported in the Theory of Wing Sections , Refer¬ 
ence [28]. Numerical values are taken from the tabulation of these experiments reported in 
Aerofoil Sections , Reference [25]. Computational geometries for the wing sections are con¬ 
structed based on governing equations and then verified to be accurate by comparison with 
measured data points in Appendix 1 of Reference [28]. The “rough” experimental results at a 
Reynolds Number of 6xl0 6 are considered. These foils are tested in the NACA two-dimensional 
low-turbulence pressure tunnel. This tunnel is specifically designed to provide an environment 
for testing relatively large wing sections in a two-dimensional environment at turbulence levels 
so low they approach the turbulence level of the atmosphere. This tunnel is further described 
in Reference [28]. 

The foils are constructed with great care to produce smooth surfaces which avoid prema¬ 
turely tripping the boundary layer near the leading edge. It is found in the testing that even 
the slightest roughness near the leading edge causes a transition to turbulence. After testing 
foils with a smooth surface, various forms of surface treatment are considered. Experiments are 
conducted using “turbulence stimulating strips” at various points along the chord. It is found 
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that their impact increases as they move forward on the foil. Experiments are also conducted 
with a sparse field of grains over the leading 8% of the chord surface. The grain density is such 
that about 5 to 10 percent of the leading edge area is coated. Grains of differing sizes are tested 
and it is shown that the initial placement of grains of any size cause the largest change then 
there is a continuing reduction in foil performance at a slower rate with increasing grain size. 
Eventually, a grain size of 0.011 inches (0.05% of the chord) is selected based on the idea that 
there is little further change in wing performance and this is still a smaller imperfection than 
would be expected by an actual wing in service. More information on the roughness selection 
process and analysis is documented in Reference [2]. 

The chord length of the foils experimentally evaluated is 2 feet. The tunnel throat is 
rectangular with a 7.5 foot height and 3 foot width. The lift coefficient is determined based on 
measured pressure distribution along the top and bottom of the tunnel. In order to evaluate 
whether or not the tunnel environment truly provided the equivalent of a two dimensional 
environment, several wing sections are selected to be evaluated at different chord lengths. Chord 
lengths of I, 2 and 4 feet are evaluated and the lift curve versus angle of attack characteristics 
are found to be consistent. In order to ensure that the zero angle of attack reference point is 
correct, lift curve versus angle of attack characteristics are measured right side up and upside 
down. 


3.1 NACA FOIL DESIGN 

3.1.1 Thickness Distribution 

The NACA 6-series thickness sections are used for this portion of the analysis. The experi¬ 
mental work demonstrates the lift curve slope to be virtually independent of the section camber 
design and highly dependent on the section thickness distribution. The design of these sections 
attempts to delay transition to a turbulent boundary layer by maintaining a favorable pressure 
gradient to a specified chord-wise location. For example, the 63012 foil section can be broken 
down as follows: 
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Figure 3-1: 6 -Series Thickness Sections 


63012 

6 = foil series designation 

3 

3 = the location of minimum pressure due to thickness is at — Chord 

0 = ideal lift coefficient is 0.0 at AOA = 0.0 
12 = maximum foil thickness is 12 % of the chord length 

The basic shape of these foils is shown in Figure 3-1. As shown in this figure, the 66 XXX 
foils have a significantly fuller aft section for the same thickness. In the comparisons which 
follow, it is observed that this change in shape results in a far more significant change in the 
experimentally measured lift coefficient than it does in the computational lift coefficient. In 
subsequent chapters of this research it is important to remember that the same experimen¬ 
tal procedure is used for both of these sections, yet the deviation between experiment and 
computations grows significantly. 
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3.1.2 Cambered Foil Design 

The camber loading of the foils in the 6-series family arises from efforts to reduce the effect of 
the “compressibility burble” on airplane wings as their speed approaches the speed of sound. 
Compressibility effects are observed more as a function of peak loading on the foil than average 
loading. The strategy used to reduce compressibility effects is to minimize the peak to average 
loading ratio by creating a constant loading over the chord of the foil. A series of foils referred 
to as the “16-series” was developed based on this premise and experimentally evaluated. The 
development and testing of these foils is reported in Reference [26]. These tests are conducted 
with an aspect ratio of six in a circular wind tunnel and extrapolated to an equivalent two 
dimensional behavior. The results show a somewhat greater over prediction than those tested 
in the NACA Langley Two Dimensional Pressure Tunnel, but have not been used for comparison 
here due to the lack of certainty of the true 2-dimensional behavior and lack of knowledge about 
whether or not the boundary layer is transitioned. 

At the time these foils were developed the understanding of the effects of viscosity on two- 
dimensional foil theory was in its infancy. The thin foil theory was applied to give a completely 
flat loading all the way from the leading edge to the trailing edge. This loading is referred 
to as the “a = 1” loading and is used as the default camber line for the extensive testing 
accomplished. It is not actually possible in a fluid with finite viscosity to maintain loading all 
the way to the trailing edge, so separation occurs and the fluid unloads the trailing edge. An 
alternative strategy used for a few of the foils was to linearly unload the foil trailing edge from a 
specified chord-wise position to the trailing edge. For example, the “a = 0.8” line has become 
a popular camber distribution for marine propellers. This camber distribution maintains a 
flat loading to the 80% chord-wise position and then linearly reduces the loading to zero at the 
trailing edge. The following is the complete designation of the 6-series family of foils: 
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66409 a = 0.8 


6 = foil series designation 

0 

6 = the location of minimum pressure due to thickness is at — Chord 
4 = ideal lift coefficient is 0.4 at AOA = 0.0 
09 — maximum foil thickness is 9% of the chord length 
a = 0.8 

foil is linearly unloaded from 80% chord position to trailing edge 
if no designation, foil design remains fully loaded to trailing edge 

Delaying of the “compressibility burble” for an aircraft wing has as a direct analogy in 
delaying cavitation in a hydrofoil. In a fluid with sufficient gas nuclei, cavitation occurs when 
the local pressure falls below the vapor pressure of the fluid in which the foil is operating. In 
the potential region of the flow around a foil, the pressure is inversely proportional to the square 
of the velocity. If the maximum velocity is minimized for a given lift coefficient, the result will 
be delaying the onset of cavitation to the maximum extent possible. 

It is based on these concepts that Brockett [6] proposed the NACA 66 modified thickness 
section coupled with the a=0.8 mean line for use in propeller design. A discussion with one 
of the engineers involved in the decision to use the a=0.8 mean line instead of the a=1.0 mean 
line indicates that the a~0.8 mean line is chosen because at the time they were not able to 
predict the behavior of the a=1.0 mean line. It turns out that the viscous effect and potential 
flow thickness effect for an a=0.8 mean line nearly cancel one another. The result is that the 
a=0.8 mean line produces lift very near the value predicted by inviscid theory. On the other 
hand, it is impossible to maintain loading to the trailing edge of the foil in a viscous fluid so the 
camber loading on a foil with a=1.0 mean line produces only about 75% of the effect predicted 
by inviscid theory. 
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3.2 IBLT Convergence Study 


MSES [10] is designed to allow analysis of multiple airfoil elements in an unrestricted domain or 
in a tunnel with walls. These capabilities make MSES particularly well suited for the analyses 
conducted in this research. 

MSES is a streamline-based Euler discretization with a two-equation integral boundary layer 
formulation. The components are coupled through displacement thickness and solved via a full 
Newton method. 

Solutions provided in this research are converged with respect to the grid geometry, grid 
density and run time parameters. Appendix A provides the details of a convergence study 
using the coupled integral boundary layer solver MSES. 

Note: The calculations here were performed with MSIS, which is an MSES version which 
assumes strictly isentropic flow outside the boundary layers. This allows considerably faster 
execution. The results were verified to be consistent with the slightly slower MSES solution. 

3.3 Parametric Modeling of Viscous Lift Prediction Errors 

In the course of analyzing the NACA foil sections it is found that there is a pattern of over¬ 
prediction of lift. The question arises whether or not something about the foil geometry or the 
flow conditions can adequately anticipate the extent to which the over prediction exists. 

In order to further explore this possibility over a thousand NACA experiments are simulated 
using MSES. A database is created that contains the experimental lift, the computed lift and 
a wide variety of boundary layer parameters near the trailing edge. This database is analyzed 
to determine whether or not the error in lift prediction can be anticipated. 

In the course of evaluating the information in this database a wide variety of parameter 
combinations are considered. It is found that in general, the size and degree of asymmetry in 
the boundary layer parameters near the trailing edge tends to result in larger lift prediction 
errors. The parameter set chosen for anticipating the error in the full viscous lift prediction is 
based on a regression analysis of the database constructed, Equation 3.1. This equation can be 
used to predict the adjustment in computational angle of attack that will be required to cause 
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Figure 3-2: NACA Foils Analysis Overview 
the full viscous computed lift to agree with the experimentally measured lift. 


[ °top ^ bottom J 

terror - ~ 

terror = Recommended change in angle of attack in radians 

8* top = Displacement thickness on upper foil surface at the trailing edge 

^bottom ~ Displacement thickness on lower foil surface at the trailing edge 

chord = Foil chord length 


Figure 3-2 shows examples of both incoherent and organized parameter selections for com¬ 
parison. On the right set of axes is the near random distribution of over prediction with a 
positive mean when compared to experimental angle of attack. On the left set of axes a plot is 
made of lift over prediction compared to the Equation 3.1 predicted lift over prediction for the 
entire database. In essence, if Equation 3.1 is superimposed onto the full viscous prediction 
the result is a significant net reduction in lift prediction errors. 

Figure 3-3 demonstrates the implications of Equation 3.1. The experimental evidence 
indicates that current viscous calculations do not adequately account for the reduction in lift 
due to viscous effects. In fact, the best fit for the data indicates that the correction is only 
about half of what it needs to be. The accuracy of this assessment and its physical cause is 
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Figure 3-3: Computed and Measured Lift Curve Slopes 

the basic foundation of the remainder of this research. 

Figure 3-4 shows a comparison of calculated versus experimental lift curve slopes for the 
family of 6-series thickness sections. This set of curves is similar to Figure 3-3, but represents 
four different thickness forms while Figure 3-3 represents only one. The lift for these foils is 
computed at several angles of attack around zero, verified to be a linear relation and the slope 
is determined based on these results. Then the angle of attack is adjusted in accordance with 
Equation 3.1 and the lift curve slope is computed again. The experimental lift curve slope is 
taken from Table 11.1 of Reference [25]. 

The top set of axes in Figure 3-4 represents the lift curve slope computed and measured 
for the 6-series foils with minimum pressure at 30% of the chord length from a thickness of 6% 
on the left to 21% on the right. The other three sets of axes represent the 6-series foils with 
minimum pressure at 40%, 50% and 60%. 

It is interesting to note that the effect of increasing thickness or moving the minim um 
pressure back on the foil has the effect of increasing the trailing edge wedge angle. It is noted 
by Hoemer [5] and other “conventional wisdom” that increasing the trailing edge wedge angle 
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can result in foil behavior which is difficult to predict. 


3.4 Conclusions 

The analyses conducted in this Chapter suggest that there is a systematic difference between 
computed results and experimental measurements. The results indicate that the differences 
between integral boundary layer theory and NACA experiments are predictable and simple. 
Some questions which are further explored in this thesis are: 

1. “What is the physical reason for these systematic over predictions in lift?” 

2. “Does the empirical fix derived based on the NACA 6-Series foils perform adequately for 
other geometries which include trailing edge thickness and other trailing edge treatments? 

3. “Does the empirical modification to the coupled integral boundary layer theory result in 
accurate local velocities as well as improved global solutions?” 
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Figure 3-4: Comparisons of Experimental and Computational Lift Curve Slopes 
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Chapter 4 


Experimental Validation 


In Chapter 3 it is shown that MSES predicts a lift force higher than that which is experimentally 
measured by NACA. In order to allow a more detailed comparison between flow codes and the 
actual flow field, measured flow velocities at strategic locations in the fluid domain are desired. 
In this chapter a set of measurements are reported and comparisons with the flow code are 
made. 

The geometries considered here have a range of features that make them of interest in 
propeller design including different forms of trailing edge modifications. They are quite different 
from the NACA geometries upon which Equation 3.1 is based, but the resulting improvement 
in flow field prediction is striking. 

4.1 Experimental Facilities 

The water tunnel at MIT is a variable pressure water tunnel described in Reference [16]. It is 
a closed-loop tunnel driven by a single 75 horsepower motor. The test section of the tunnel is 
twenty inches square and approximately four feet long. The maximum achievable test section 
velocity is 30 Upstream of the test section is a five to one contraction fitted with fine 

honeycomb mesh and a wake screen to promote flow uniformity. 

The test section has removable plexiglass walls on all four sides for ease of access and to allow 
fluid velocity measurements via laser Doppler velocimetry. For this experimental technique 
the velocity is inferred based on the Doppler shift in light frequency scattered by a neutrally 
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buoyant seed particle in the fluid. Detailed error analyses are conducted in Reference [21]. 

The flow conditions at the mid-span of the test section is demonstrated to be a two dimen¬ 
sional flow environment in the following ways: 

1. In the NACA low turbulence wind tunnel a test section with a width of 3 feet is used 
to perform experiments on foils with a chord length of two feet. These experiments are 
shown to provide equivalent two dimensional performance at the mid-span based on the 
fact that the lift versus angle of attack curves are consistent for experiments conducted 
on foils with a chord of six inches, one foot and two feet [28]. The ratio of the foil chord 
length to tunnel span is similar for the NACA experiments and those conducted at the 
MIT variable pressure water tunnel. 

2. In a thesis by Jergen Jorde [15], fluid velocities are measured across the span of the tunnel 
and are shown to be constant outside of the front and back side boundary layers which 
are less than one inch thick. It is demonstrated that the circulation as a function of span 
is constant outside of the tunnel wall boundary layers. 

3. In a thesis by Randall Fairman [11] it is demonstrated that boundary layer migration 
from the tunnel walls does not effect the flow field at the center of the tunnel. This 
is accomplished using a number of fence designs to prevent three dimensional boundary 
layer migration. 

4.2 Experimental Procedures 

Seven hydrofoil geometries are tested. The hydrofoils tested range from fourteen inches to 
eighteen inches in chord and span the entire width of the tunnel. The Reynolds Number for 
these foils are 1.4 x 10 6 to 2.0 x 10 6 respectively. The boundary layer on the foils is tripped 
at the chord position. The development of these foils and turbulence strips is found in 
Reference [24] and Reference [18]. 

Each foil is mounted at the midpoint of the test section on a shaft which allows the foil to 
rotate about its axis. Each foil is also supported near the trailing edge to provide rigidity. For 
each geometry given, the axis of rotation is at (0,0) and the axis located vertically in the center 
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of the tunnel. The foil geometries in the Appendices are given at a zero angle of attack. 

The actual position of the surface of the foil is measured via laser crossing at each ex¬ 
perimental angle of attack. In general, the laser positioning device is capable of resolving 
position within 0.01 mm. The foil surface is rotated in the computational domain until the 
computational and experimental non-dimensional foil surface positions agree. This procedure 
for measuring angle of attack is found to be repeatable within 0.05 degrees. Lift curve versus 
angle of attack is measured both right side up and upside down for the “HRA Blunt” foil. 
These measurements demonstrate that the error in angle of attack is less than 0.1 degrees. 

The tunnel velocity is held constant within 0.1% of free stream velocity. With a sample 
size of 2,500 valid velocity data points at each spatial location it is found that the velocity 
measurements are repeatable to within 0.01 The flow is non-dimensionalized by the 

measured inflow velocity for the given experiment which is typically about 5^. The resulting 
error band for measured pressure coefficient is ±0.004 (Equation 4.1). 


C p 

C p 

Ulocal 


Vlocal 


u 0 


0.5 


^local + V Uxxd 

2 uz, 


where 

pressure coefficient 

horizontal velocity measured at x,y 

vertical velocity measured at x,y 

average horizontal velocity at inlet side of control volume 


(4.1) 


4.3 HRA Series of Foils 

The HRA series of foils are designed in conjunction with research seeking to improve cavitation 
performance in marine propellers. They are designed to obtain a flat pressure distribution 
over the majority of the foil with reduced loading near the leading and trailing edges. The 
reduced loading near the leading edge is intended to allow improved cavitation performance at 
off design angle of attack conditions. The details of the design and analysis of these foils can 
be found in Reference [18]. 
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4*3.1 HRA Baseline Foil 



Figure 4-1: HRA Fine Trailing Edge Foil Data 

The HRA baseline foil has a fine streamlined trailing edge. The foil geometry and tabulated 
experimental results are found in Appendix B. The foil has a chord length of 16 inches which 
yields a non-dimensional tunnel height of 1.25 chord lengths. 

The fluid domain is analyzed using MSES as shipped and with the change in infinite angle 
of attack implied by Equation 3.1. The tunnel walls are modeled as a slip surface. In a subse¬ 
quent section it is shown using the RANS program CFDSHIP that the difference in calculated 
velocities between a slip tunnel wall and a no-slip tunnel wall is insignificant compared to vari¬ 
ations under consideration here. Results are plotted as “Experimental”, “MSES calculated” 
and “MSES adjusted” respectively. In this case the foil drag is relatively low and the resulting 
displacement thicknesses are relatively small. 

Some noteworthy observations from Figure 4-1: 

1. At all three experimental angles of attack the calculations as modified by Equation 3.1 
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provide superior results. 

2. The implied error in lift coefficient at the highest angle of attack is ~5% for the “MSES 
calculated” results and less than 1% for the “MSES adjusted” results. 

3. The adjustment in infinite angle of attack improves not only the global force prediction, 
but also the local velocity predictions. 

4.3.2 HRA Blunt Foil 





Figure 4-2: HRA Blunt Trailing Edge Foil 

The HRA Blunt Foil is designed to have the same pressure profile as the HRA fine trailing 
edge foil but features a 2% thick trailing edge. The foil geometry and tabulated experimental 
results are found in Appendix C. The foil has a chord length of 16 inches which yields a 
non-dimensional tunnel height of 1.25 chord lengths. Figure 4-2 shows the computations and 
measurements for this foil. 

Observations: 

1. The presence of a large trailing edge thickness does not effect the magnitude of the over 
prediction. 
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2. As the trailing edge boundary layer approaches the symmetric condition the raw calcula¬ 
tion, adjusted calculation and experimental results converge. 

4.3.3 HRA Beveled Foil 





Figure 4-3: HRA Beveled Trailing Edge Foil 

The HRA. beveled trailing edge foil is identical to the HRA blunt trailing edge foil with the 
exception of a beveled suction side over the last 2.5%. This feature is shown to cause the vortex 
shedding pattern to become significantly asymmetric in Figure 1-1. Figure 4-3 is a summary 
of the data for this foil and the results are presented in tabular form in Appendix D. 

Observations: 

1. The presence of the beveled section near the trailing edge results in virtually no change 
in experimentally measured lift, yet yields a significant increase in the predicted lift. 

2. The change in infinite angle of attack suggested by Equation 3.1 results in a significant 
improvement in prediction performance. 

3. Even after adjusting the infinite angle of attack there remains a slight over-prediction in 
lift coefficient primarily caused by inaccuracies near the trailing edge. 
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4. In the MSES adjusted solution there is a slightly lower than actual loading near the 
leading edge and a slightly higher than actual loading near the trailing edge. When the 
MSES adjustment in accordance with Equation 3.1 is in error, this seems to be the nature 
of the error. It appears this error is amplified by flow curvature near the trailing edge. 
In each case here it seems that the overall lift force is still well predicted by Equation 3.1. 
The local loading could be more adequately predicted by coupling a more slight change 
in angle of attack with a modification to the “Kutta Condition” near the trailing edge. 
It is suspected that asymmetric vortex shedding results in pressure biasing which cannot 
be captured by an isotropic turbulence model. These issues will be further addressed in 
subsequent sections. 

5. At the largest angle of attack the implied error in lift coefficient is about 0.16 (« 20%) 
for the basic MSES calculation. 

4.4 B1 Series of Foils 

The B1 series of foils are constructed based on the Brockett [6] modification to the NACA 66 
thickness section. These foils are commonly used in propellers so the experimental character¬ 
istics are of particular interest to the Naval Architect. The details of the design and analysis 
of these foils can be found in Reference [24]. 

4.4.1 B1 Blunt Foil 

The B1 Blunt foil has a modestly blunted trailing edge (~ 0.75%) used for analysis as an 
alternative to the standard anti-singing design. Figure T4 is a summary of the data presented 
for this foil in tabular form in Appendix E along with the MSES calculations. 

Observations 

1. As in the case of the HRA Blunt foil, the performance of Equation 3.1 in predicting the 
computational error is good. 

2. The implied error in lift coefficient for the basic MSES calculations is ~10%. 
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Figure 4-4: Bl Blunt Trailing Edge Foil 

4.4.2 Bl Anti-singing Foil 

This foil has a trailing edge which is modified in accordance with industry best practices used 
to prevent propeller “singing”. Figure 4-5 is a summary of the data presented for this foil in 
tabular form in Appendix F along with the MSES calculations. 

Observations 

1. This foil has slightly greater flow curvature near the trailing edge than the Bl Blunt foil 
and less than the HRA Beveled trailing edge foil. It also has a slightly greater over 
prediction of differential pressure at the trailing edge than the Bl Blunt foil and less over 
prediction of trailing edge differential pressure than the HRA Blunt foil. 

2. At times the similarity in profile shapes between the computations and the measurements 
can cause one to doubt the importance of the differences. Consider: 

(a) This foil is designed in accordance with typical best practices for military propeller 
blade shapes. 

(b) The experiment operating at an angle of 2.25 degrees is well within the expected 
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Figure 4-5: B1 Antisinging Trailing Edge Foil 
operating region. 

(c) The magnitude of the over prediction in lift at this angle is over 10%. 

4.4.3 B1 Cupped Foil 

This foil has a radically cupped trailing edge used to evaluate the capabilities of computational 
fluid dynamics in capturing the effects of trailing edge details on the global flow field. Figure 
4-6 is a summary of the data presented for this foil in tabular form in Appendix G along with 
MSES calculations. 

Observations 

1. This foil was originally designed to highlight the weaknesses of computational fluid dy¬ 
namics. The difficulties in evaluating the mean lift on this foil drove the investigation 
into other experimental efforts throughout the world and at MIT. The investigations of 
other experiments and computations led to Equation 3.1 which is now applied back to 
this foil. The result is an accurate lift prediction for this foil. 

2. In terms of velocity distributions, this foil has greater trailing edge curvature and a greater 
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Figure 4-6: B1 Cupped Trailing Edge Foil 

over prediction of pressure loading near the trailing edge when compared to the B1 Blunt 
foil. 

4.5 Conclusions 

4.5.1 It is demonstrated: 

1. A modification to the infinite angle of attack suggested by Equation 3.1 provides not 
only improved predictions for the global lift coefficient, but also an improvement in local 
velocity predictions. 

2. The presence of trailing edge thickness does not effect the ability of MSES to predict lift 
coefficient. 

3. Equation 3.1 is equally effective at improving MSES predictions with or without trailing 
edge thickness. 

4. In the presence of high trailing edge curvature there remains an over prediction in trail¬ 
ing edge loading and an under-prediction in leading edge loading that results in a good 
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estimation of global lift, but results in errors in predicting the moment coefficient for a 
wing section. 

4.5.2 Questions for further exploration: 

1. What is the physical cause of these errors? 

2. Can a different form of turbulence modeling such as RANS produce superior results? 
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Chapter 5 


Integral Boundary Layer Theory 
Compared to RANS 


The coupled integral boundary layer code MSES is specifically designed for two dimensional 
fluid dynamics analyses associated with foils. The process of becoming an adequate user of this 
flow code required an investment of about a month to really understand the inputs, outputs, 
and steps necessary to verify convergence. After this initial investment a full convergence study 
for a given case from grid generation to a final solution takes about thirty minutes. 

For this research the primary RANS code utilized is CFDSHIP. CFDSHIP is a fully implicit 
RANS code developed with the stated goal of modular code and ease of user modification for 
research. A full report on the details of the code development can be found in Reference [27]. 
I have personally used or attempted to use IFLOW, DTNS[13], UNCLE [3] and CFDSHIP. In 
my opinion, the CFDSHIP development team has accomplished the goal of code simplicity and 
has given the user the capability to freely add subroutines and modify boundary conditions 
without excessive danger of unknown effects. Within three months of starting to use the code 
the various modules and subroutines became intuitive. I spent significant time with each of 
the other codes as well but did not find the codes to be as organized or understandable. 

The process of becoming an adequate RANS user is, for this author, at least a two year 
process. The complications of generating an adequate grid with high density in the appropriate 
locations and then determining run time parameters that will allow a RANS code to run to 
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convergence has been a tedious process at best. In the end, the result is a solution which is 
independent of grid geometry and topology. Having completed the initial investment, it now 
requires approximately two days to bring a single foil to convergence and verify that the solution 
is independent of grid density and topology. 

As a note of interest, I did find that there were certain grid topologies combined with 
particular grid densities that would provide a solution that matched the experiment quite well. 
On each of these occasions careful grid refinement with proper convergance resulted in solutions 

which were again essentially identical to the MSES solutions. 

It should be noted that although RANS is far more difficult to use, it is also far more 
general in its ability to handle arbitrary geometries. MSES is specifically designed to analyze 
a two dimensional foil like geometry. CFDSHIP is designed to be effective in a complex three 
dimensional environment associated with analyses of the types of flow that occur in a marine 
environment. The increased time and complexity associated with RANS described here is in 
part necessary because of the increased generality of the flow solutions. 

5.1 Grid Topology 

Figure 5-1 shows the overall grid geometry used and the leading edge details. A single zone 
was chosen for the entire inflow region in order to minimize the number of zone boundaries 
required. It was found that the presence of zone boundaries significantly increased the time 
to convergence and the potential for instability while marching toward a solution. With this 
geometry the tunnel walls and foil leading edge are all able to support a high grid density while 

using a structured grid. 

The zones above and below the foil have increased density near the leading and trailing 
edges where the flow velocities have the highest gradients. These zones extend into the leading 
and trailing edge zones in order to maintain higher order accuracy across the zone boundaries. 
Experience indicates that if the upper and lower zones extend into the leading and trailing edge 
zones with an overlap of ten cells there is a significant improvement in stability and there is 

continuity of higher order terms across the boundary. 

Figure 5-2 shows the details of the trailing edge grid topology. Here a single zone is again 
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Figure 5-1: Leading Edge Grid Topology 

used to enhance stability. A small aspect ratio (ten or less) is used at the trailing edge in order 
to ensure that the large flow gradients in this region can be captured. 

5.2 Boundary Conditions 

The tunnel walls can be modeled as slip or no-slip surfaces. MSES models the walls strictly 
as slip surfaces so the comparisons with MSES use slip surfaces. The experiment has a real 
wall so a no slip surface analysis is also required. It is found that the local flow conditions 
around the foil remain mostly unaffected by the wall boundary condition chosen. The pressure 
gradients at the walls are set to zero. 

The inflow condition applies a horizontal velocity of unity and a zero gradient in the hori¬ 
zontal direction for vertical velocity and pressure. Other inflow possibilities such as constant 
pressure and zero vertical velocity were used for further evaluation. These other inflow condi¬ 
tions provide results consistent with those presented here. 

The exit boundary condition is modeled as follows: 
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Figure 5-2: Trailing Edge Grid Topology 


1. The pressure gradient in the horizontal direction is set to zero. 

2. The velocity gradients in the horizontal direction are set to zero. 

3. Mass flow out of the flow domain is calculated. 

4. A “differential” horizontal velocity is uniformly applied to the outflow boundary in order 
to force global conservation of mass. 

It has been our experience at the MIT water tunnel that RANS calculations for internal flow 
problems are quite often extremely unstable. There have been many attempts at exit boundary 
conditions. The set of exit boundary conditions used here have been applied to internal 
flow calculations including waterjets, propellers and foils with overall stable results. The 
RANS codes CFDSHIP and DTNS have both been exercised with these boundary conditions. 
The effect of the differential horizontal velocity is to maintain a reasonable solution during 
the convergence process and ultimately approach a zero differential horizontal velocity when 
convergence is achieved. 







Iteration Number 

Figure 5-3: Iterative Convergence History 

As mentioned above, pressure gradients are used as boundary conditions on all surfaces. 
This makes it necessary to fix the pressure at one location in the fluid domain. For these 
analyses the pressure at about the center of the upstream domain is set to zero. 

5.3 Iterative Convergence 

It is imperative to ensure that each case is fully converged with respect to the number of iter¬ 
ations. The most reliable method of verifying iterative convergence is to track the “residuals” 
and a relevant physical parameter. The residuals are expected to continue to decrease and 
the relevant physical parameter should become steady. For this study the fluid dynamic lift 
was measured each iteration by applying a numerical contour integration to a path in the fluid 
domain that surrounds the foil. When the predicted lift becomes constant and the residuals 
continue to decrease the problem is considered to be adequately converged. Figure 5-3 shows 
the convergence history for one of the B1 Blunt foil cases. 
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Figure 5-4: Y Plus Convergence for B1 Blunt Foil 

5.4 Grid Density 

5.4.1 Y Plus Convergence 

The first issue explored is the height of the first grid node from the surface of the foil. The 
non-dimensional height above the foil surface is commonly referred to as the “Y + value.” Most 
published results indicate that a Y + value of one or less is needed to ensure the details of the 
boundary layer are captured. Equation 5.1 is the basis for calculating this parameter. 


Y+ = — (5-1) 

v 

y = distance to foil surface 
v = kinematic viscosity 

^ M 

t w = wall shear stress 
p = density 
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Figure 5-5: Bulk Grid Area Maximum Spacing 

One of the practical problems in selecting a grid is that the wall shear stress is not known 
until after the solution is known. The experience of this process is that if the first node is 
too large it will give an inaccurate solution and if the first node is to close the solution will 
become unstable. The manner in which this dilemma is solved is to start off with a first node 
Y + value slightly too large and then reduce the first node height until no further changes in 
the solution occur. Figure 5-4 shows the lift coefficient changes as the Y+ value is reduced to 
ensure convergence. It should be noted that each of the CFDSHIP Predicted values requires 
about two hours of computer time on a 1.8 GHz AMD Athlon Processor. 

5.4.2 Grid Spacing 

The upper and lower zones are separated into regions according to anticipated flow gradients. 
The leading edge region of the foil extends upstream and downstream of the leading edge by 5% 
of the chord length. The trailing edge region extends extends 5% upstream and downstream 
of the trailing edge region. The foil region extends 10% above and below the foil. 

Figure 5-5 shows the lift coefficient verses maximum grid dimension in the bulk of the fluid. 
In each case this needs to be verified, but in cases analyzed for this research a bulk grid spacing 
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Figure 5-6: B1 Cupped Foil Geometry 


of 4% of the chord length is adequate. A similar process must be conducted to verify that the 
grid spacing in the regions of higher flow gradients is also converged. The degree of refinement 
required in these regions varies depending on the details of the flow. Typical grid spacing in 
these regions is between 0.05% of the chord length for the B1 Cuppped foil and 2% of the chord 
length for the HRA Fine Trailing Edge Foil. 

5.4.3 Turbulence Models 

CFDSHIP is programmed to support the K-e, K-w and Baldwin/Lomax turbulence models. 
DTNS is programmed to support the K-e and Baldwin/Lomax turbulence models. The foils 
under consideration here have been evaluated using the CFDSHIP K-w turbulence model, the 
DTNS K-c turbulence model, the DTNS Baldwin/Lomax turbulence model and the MSES 
integral boundary layer turbulence model. In each of these cases it is found that the details 
of the local boundary layers may vary, but the global lift coefficient as indicated by the flow 
velocities in the potential flow region where measurements were taken is almost identical. 


5.5 Foil Comparisons 

In Chapter 1 the results of a previous blind test on the Bl-Cup foil were discussed as a part of 
the motivation for this research. In this blind test several authors of several codes (including 
those listed here) evaluated the B1 Cupped foil shown in Figure 5-6. The results reported 
were so varied that the conclusion was drawn that the details of the turbulence modeling or 
the software were inconsistent. The conclusion of this research is that the most likely reason 
for the previous variation in predictions is inadequate grid convergence study prior to drawing 


conclusions. 




X/C Position 

Figure 5-7: Solution Comparisons for the B1 Cupped Trailing Edge Foil 

Figure 5-7 shows a fully converged solution from the CFDSHIP K-lj model, a fully con¬ 
verged solution from MSES and the experimentally measured datapoints. This foil experiences 
significant vortex shedding, significant separation prior to the trailing edge and has a very thick 
boundary layer. Shown here are the two most radically different turbulence models under 
consideration: IBLT and a two equation eddy viscosity model. The similarity in the results is 
striking. 

The convergence procedure and comparisons similar to Figure 5-7 are conducted and the 
consistent result is that IBLT and RANS provide virtually identical results. The codes have 
been verified to give the same result for the HRA fine trailing edge foil, the HRA blunt trailing 
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edge foil, the B1 Blunt trailing edge foil and the B1 cupped trailing edge foil. In every case 
considered there is an over prediction of lift coefficient approximately consistent with Equation 
3.1. 

5.6 RANS Solution with Walls 

One consideration is that the growth of a boundary layer on the upper and lower walls cause 
an overall acceleration of the fluid in the tunnel. In order to evaluate the magnitude of this 
concern, the water tunnel is modeled as a no-slip surface with a boundary layer. It is found 
that there is a slight increase in the velocity from the leading edge to the trailing edge. This is 
shown in Figure 5-8 by a slightly lower pressure coefficient in the back part of the foil. There 
is virtually no change in the overall pressure loading and it is certain that this phenomenon is 
not the source of the prediction errors with either RANS or IBLT. 

5.7 Conclusions 

This chapter demonstrates that RANS and coupled IBLT give essentially identical lift coefficient 
and pressure profile predictions when fully converged. With this established the following 
conclusions can be drawn: 

1. These codes are developed completely independently, so it is certain that the differences 
between experiments and predictions are not caused by errors in the implementation of 
the underlying assumptions about the physics. 

2. The differences between the codes and experiments can be predicted based on flow condi¬ 
tions which could be expected to effect the accuracy of the computational domain, not the 
experimental results. The differences between the codes and experiments are consistent 
across a spectrum of experimental facilities and geometries which make them unlikely to 
be caused by experimental error. 

3. The codes use completely different methods to model the underlying assumptions about 
the physics and yet reach the same conclusion; thus, it is concluded that there is likely to 
be a common assumption which is in error. 
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The investigation into the underlying assumptions about the physics being modeled is ex¬ 
plored in subsequent chapters. 
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Chapter 6 


Analysis of Assumptions 


In Chapter 5 the conclusion is drawn that the most probable explanation for the observed over 
prediction in lift is an inadequate assumption about the underlying physics which is common 
to both IBLT and steady RANS codes. In Chapter 2 an overview of these two approaches to 
solving the fluid dynamic problem is discussed. It is appropriate here to further consider the 
derivation of the underlying equations in order to highlight potential erroneous assumptions. 


6.1 Navier-Stokes Equations 

The Reynolds Averaged Navier Stokes (RANS) Equations are a time average of the Navier- 
Stokes Equations; therefore, the assumptions used in deriving the Navier-Stokes Equations are 
inherent in the RANS formulation. The Navier-Stokes Equations are derived as follows: 

First, Newton’s second law (F — ma ) is applied to a fluid particle. Then both sides are 
divided by the volume and rearranged in accordance with tradition. This gives Equation 6.1. 
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f — fbodyg “t" fbodymhd ~b fsurface (®*^) 

density 

particle derivative in Eulerian system 
gravity forces 

magneto hydrodynamic forces 
surface forces 

Assuming the magneto hydrodynamic forces are not relevant and applying a standard stress 
tensor notation gives the common form shown in Equation 6.2. 

DV (a 

p— = pg + V • Tij (6.2) 

To this point in the derivation, the only assumption made is that the magneto hydrodynamic 
forces are not relevant, which seems safe. At this point in the derivation, assumptions about 
the fluid properties are required. Three basic assumptions about the fluid properties were 
made to obtain Stokes Hypothesis: 

1. The fluid is continuous and the stress tensor (Tij) is a linear function of the strain rate. 

2. The fluid is isotropic. 

3. When strain rates are zero, deformation laws must reduce to the hydrostatic condition. 

These assumptions have been shown to be accurate for all gases and most common fluids, 
including water. In this case the fluid under consideration is water, so the assumption of an 
incompressible fluid is made. The direct result of stokes hypothesis and an incompressible fluid 
is Equation 6.3 which gives the deformation law for incompressible Newtonian fluids. 


DV 

p 

p 

DV 

Dt 

fbodyg 
fbodymhd 
fsurface 


'Tij — P^ij M 


/ dui duj 
\ dxj dxi 


) 


(6.3) 
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Combining Equation 6.1 and 6.3 gives the incompressible Navier Stokes Equations. Equa¬ 
tion 6.4 is shown for the x direction instead of using indicial notation for clarity. 


Du dp d 

P ~Dt ~ P9x ~dx + dx 


( du\ d 

l 2 ^) + Ty 


{du dv\ d f dw du 
p \dy^dx) dz P \ dx + dz 


6.2 Reynolds Averaged Navier-Stokes Equations 


The Navier-Stokes equations are a direct result of Newton’s Second Law and therefore apply 
equally in steady state or in a transient. To the extent that the assumptions made about the 
fluid properties are accurate, the Equations are accurate at any instant in time. 

The concept employed by Reynolds is to model the fluid velocities and pressures with an 
average term and a fluctuating term. There are no further assumptions about the fluid or 
the physics. The result of time averaging is most easily seen by replacing the deformation 
law which was a result of stokes hypothesis about the fluid properties (Equation 6.3) with a 
deformation law that includes the effect of unsteadiness, Equation 6.5. 

r,i = -pSij + » (§^ + fg) - (6.5) 

p ^ 7 ” + ^ = steady (laminar) part of the stress 

pu'ju'j = unsteady (turbulent) part of the stress 


It is important to note that the time averaging of the Navier Stokes Equations has added 
no new assumptions. To this point, the only assumptions made concern the fluid properties 
in accordance with Stokes Hypothesis. These assumptions seem unlikely to be the source of 
observed errors. 

While it is true that no new assumptions have been made, it is also true that nine new 
variables have been introduced. These nine new variables bring with them no real physical laws 
such that a great need for empirical modeling results. There have been extensive discussions 
among experts in the field concerning the possibility that these nine new variables sufficiently 
increase the difficulty in solving the physics that perhaps a totally new direction should be 
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taken. The following quote is from White, a noted expert in fluid dynamics: (Reference [29]): 


Although statistical theory and numerical simulation are viable options, most of the 
research on turbulent-flow analysis in the past century has used the concept of time 
averaging. Applying time averaging to the basic equations of motion yields the 
Reynolds equations, which involve both mean and fluctuating quantities. One then 
attempts to model the fluctuation terms by relating them to mean properties or their 
gradients. This approach may now be yielding diminishing returns: Lumley (1989) 
gives a stimulating discussion of how time averaging might outlive its usefulness. 

The Reynolds equations are far from obsolete, however, and form the basis of most 
engineering analyses of turbulent flow. 

6.3 Turbulence Modeling 

The concept of turbulence modeling is, in its very essence, simply making an attempt to deter¬ 
mine the nine unsteady terms in Equation 6.5. The difficulty lies in the fact that there are no 
further physical laws available. The following quote is from Reference [29] and highlights the 
difficulty associated with these terms: 

Thus the mean momentum equation is complicated by a new term involving the tur¬ 
bulent inertia tensor u',u',. This new term is never negligible in any turbulent flow 
and is the source of our analytic difficulties, because its analytic form is not known a 
priori. In essence, the time-averaging procedure has introduced nine new variables 
(the tensor components) which can be defined only through unavailable knowledge 
of the detailed turbulent structure. The components of u^Uj are related not only 
to fluid physical properties but also to local flow conditions (velocity, geometry, sur¬ 
face roughness, and upstream history), and no further physical laws are available to 
resolve this dilemma. In a two dimensional turbulent boundary layer, the only sig¬ 
nificant term reduces to u'v f , but even this single term requires extensive scratching 
about to achieve an analytic correlation which is semiempirical at best. Some of 
the empirical approaches have been quite successful, though rather thinly formulated 
from nonrigorous postulates. 
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Earlier in this chapter it is concluded that the RANS Equations in and of themselves contain 
few assumptions that can be questioned. The fact that the RANS and IBLT approaches provide 
virtually identical results gives the conclusion that it is very unlikely that the problem lies in 
the implementation. Given these observations it seems likely that source of the observed over¬ 
predictions must be in the turbulence modeling. This conclusion is buoyed by the fact that 
the quantity which best predicts the magnitude of the over-prediction is in its very essence a 
measure of the magnitude and assymetry of the local turbulence at the trailing edge of the foil. 
On the other hand, this research has implemented several different turbulence models that all 
provide consistent solutions. 

6.4 Conclusions 

There are a long list of experiments at multiple facilities that all demonstrate that lift is over¬ 
predicted by the flow codes. On the other hand, there are four different flow codes which use 
two radically different methods to model the fluid physics. There are four different methods 
used to model the turbulence. All of these flow solutions provide the same answer which is 
consistently higher than the experiments by a predictable amount. The dilemma forces us to 
stand back and ask the question: 

Is there something about the flow around a foil which is unmodelled? 

The answer to this question is yes! 

1. All of the flow codes considered calculate only the time averaged behavior of the flow 
around the foil. The reality is that each time a discrete vortex is shed from the foil into 
the wake there is a surge in global lift coefficient. This surge in global lift coefficient 
results in another vortex of the opposite sign being shed and the process repeats. In a 
time averaged code, there are no discrete vortices shed and therefore there is no global 
response. The manner in which this phenomenon could manifest itself in a lift over¬ 
prediction would be through errors in wake trajectory caused by turbulence parameter 
inaccuracies. 
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Figure 6-1: Cupped Foil Near Wake Turbulence Characteristics 

2. All of the flow codes considered do not account for non-isotropic turbulence characteristics. 
Measurements using particle image velocimetry indicate that the flow around the trailing 
edge of the foils with large over prediction of lift experience highly non-isotropic turbulence 
characteristics. Figure 6-1 shows the measured horizontal and longitudinal turbulence 
characteristics for the B1 Cupped Foil at an angle of attack of about 0.5 degrees. Notice 
that at one percent of a chord length downstream of the foil the horizontally oriented 
turbulence intensity is about four times higher than the vertically oriented turbulence 
intensity! 

3. It is likely that the former observation (unsteadiness) at least in part results in the latter 
observation (non-isotropic turbulence). In any case, both of these observations could 
explain how the flow codes could develop the same solution and still be different from the 
experiments. Other researchers have run into trouble predicting the characteristics of 
complicated flows like this. The following quote is from Henkes [14]: 

An example of a complex flow is turbulent boundary-layer separation. The sep¬ 
aration shows strong curvature and has strong velocity gradients. Such complex 
flow types can more accurately be predicted by differential Reynolds-stress mod¬ 
els, which solve a partial differential equation for each component of the Reynolds- 
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stress tensor. For example, to calculate a 3D flow 6 equations are solved, (e.g. 
u a ,v' 2 ,w' 2 ,u'v',v'w' and v'w', respectively). In addition, the model also solves 
a differential equation for the dissipation rate. The reason that the differential 
Reynolds-stress models are most suited for complex flows is that the production of 
turbulence energy P does not have to be modelled. Therefore the production due to 
simple shear, curvature, 3D effects and rotation is represented exactly. 

For flows that contain relatively large-scale spatial structures and low-frequency un¬ 
steadiness, line near wake flows, flow buffeting and dynamic stall, the use of a 
Large-Eddy Simulation seems to be a promising approach. In contrast to the ear¬ 
lier mentioned models, which all try to model the time-averaged equations, the LES 
only models the small-scale structures, whereas the large-scale structures are calcu¬ 
lated from an unsteady time integration. As the small-scale turbulence is almost 
isotropic, the so-called subgrid models in LES are expected to be not very crucial. 
LES was originally developed for meteorological applications, like the earth’s atmo¬ 
spheric boundary layer, which typically consists of large-scale structures. More 
recent research on LES is also devoted to aeronautical boundary layers. 
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Chapter 7 


Unsteady Analysis 


In the previous chapter it is concluded that the most likely improper assumption used in the 
steady flow codes considered is related to the fact that discrete vortices are shed at a discrete 
frequency related to the Strouhal Number at the foil trailing edge. These vortices yield a 
non-isotropic turbulence with a relatively low frequency variation that cannot be anticipated 
by a typical isotropic turbulence model. 

7.1 Observed Vortex Shedding 

In the introduction, reference was made to high speed photography which observed vortex 
shedding behind foils with blunt and beveled trailing edges (Reference [18]). These foils 
exhibited a highly organized vortex shedding pattern shown in Figure 7-1. The total trailing 
edge thickness in the case of the blunt foil is two percent, so the wavelength of the pattern 
observed is on the order of ten percent. 

An experimental method called particle image velocimetry can be used to obtain an instan¬ 
taneous velocity profile over a region of the fluid domain. This experimental method uses two 
pictures of the flow field very close together in time and then uses the distance the particles 
travel and the difference in time between pictures to obtain velocities. This experimental 
method is applied for the B1 Cupped foil (Section 4.4.3) and the HRA fine trailing edge foil 
(Section 4.3.1) in order to evaluate the characteristics of the vortex shedding. 
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Figure 7-2: Wake Observation for the HRA Fine Trailing Edge Foil 

7.1.1 HRA Fine Trailing Edge Foil Vortex Shedding 

Figure 7-2 shows a single image of the wake from the particle image velocimetry. The contours 
are based on the magnitude of the horizontal velocity and the streamlines are integrated from 
the raw horizontal and vertical velocities. The estimate of wavelength indicated is based 
primarily on the contours of fluid velocity. There are fluctuations in the angle of the trailing 
edge streamlines of about ±5 degrees depending on the phase at which the pictures are taken. 
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Figure 7-3: Wake Observation for the B1 Cupped Trailing Edge Foil 

7.1.2 B1 Cupped Trailing Edge Foil Vortex Shedding 

Figure 7-3 shows a single image of the wake of the cupped foil from the particle image velocime- 
try. The contours are based on the magnitude of the horizontal velocity and the streamlines 
are integrated from the raw horizontal and vertical velocities. It is noted that there is a signifi¬ 
cantly larger wavelength than for the HRA fine trailing edge foil and that the vortices continue 
impact the streamline direction much further downstream. There are very large fluctuations 
in the angle of the trailing edge streamlines. Based on observing several hundred pictures such 
as this it is noted that the pressure side vortices (counterclockwise revolution) have a strong 
tendency to have a very small core and a large local effect on velocity. The suction side vortices 
(clockwise revolution) have a tendency to have a much larger core and a much less distinct effect 
local velocities. 

Kelvin’s theorem indicates that the time averaged vorticity shed from the foil must be zero 
(neglecting the starting vortex). Figure 6-1 shows the result of time averaging 2,500 of the 
instantaneous particle image velocimetry data fields. In this Figure there is a strong anisotropic 
nature to the turbulence on the pressure side of the foil and somewhat less on the suction side. 
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It seems that the most likely explanation for the difference from the suction to pressure side 
anisotropic turbulence characteristics has to do with the effective core size of the vortices being 
shed. On the suction side of the foil there is a much larger boundary layer and greater formation 
time starting somewhat upstream of the trailing edge such that when the vortex leaves the foil 
it has a relatively large core. On the other hand, the pressure side boundary layer is very thin 
and the vortex is actually shown to form just downstream of the trailing edge of the foil and is 
very “tightly wound”. The result is a highly asymmetric vortex shedding behavior. 


7.2 Asymmetric Vortex Shedding Analysis 

The particle image velocimetry observations make it appear that there is a clear distinction 
between the local turbulence and a larger low frequency turbulence associated with the vortex 
shedding. It is hoped that by gaining insight into the vortex shedding effects the problem can 
be better understood. 
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(7.1) 


Equation 7.1 is derived in Reference [29] and gives the velocity profile of a vortex which 
starts as a point vortex at time zero and then has a core growth and viscous dissipation over 
time. A small Fortran program is written to evaluate how a set of discreet vortices being shed 
from the trailing edge of a foil will effect the flow field. This program is provided in Appendix 
H. 

First of all, it is interesting to see the essence of how the RANS Equations attempt to model 
a vortex field. Figure 7-4 shows a set of symmetric vortices in a typical time accurate format 
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Figure 7-4: Time Averaged vs Time Accurate Representations of Symmetric Vortices 


that you would see in a picture and also in a time averaged format. In the case shown in Figure 
7-4 the vortices are symmetric and counter rotating such that the net effect on the outflow angle 
is zero. The vortices each have a vortex strength equal to one tenth of that of the B1 Cupped 
foil and the wavelength is based on the observed wavelength from PIV for the B1 Cupped foil. 

The next question considered is how an asymmetry in core size effects the mean flow field. 
The reference starting time of the vortices is adjusted and then the calculations are performed 
again. The reference starting time is used to give the vortex an initial core size similar to the 
boundary layer thickness on the associated side of the foil. Figure 7-5 shows the time accurate 
and time averaged effects of asymmetry in the vortex core size. In this simulation the vortices 
are started at the same horizontal location, but with different reference formation times. It is 
interesting to note that the asymmetry in the RMS velocities shows the same character as the 
asymmetry in the experiment (Figure 6-1). It appears that the issue of vortex core size may 
account for the observed turbulence parameter asymmetry. 

Now the question arises how the flow field would be effected if the vortices were released 
from slightly different horizontal positions. For example, on the cupped foil it appears that the 
pressure side vortex forms just barely beyond the trailing edge of the foil and the suction side 
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Figure 7-5: Time Accurate and Time Averaged Effects of Core Size Asymmetry 

vortex forms about one to three percent upstream of the trailing edge. Figure 7-6 shows the 
time accurate and time averaged effect of releasing the pressure side vortex at the trailing edge 
and the suction side vortex two percent upstream. With this type of asymmetry the suction 
side vortex has a greater time that it effects the flow field and there is an obvious net change 
in the angle of attack. This type asymmetry explains why an adjustment in the foil angle of 
attack in accordance with Equation 3.1 results in improved velocity predictions throughout the 
flow domain. 

In summary, these analyses indicate that asymmetric vortex shedding can effect the outflow 
angle and cause non-isotropic turbulence. These observations are consistent with experimental 
observations and could explain the offset in angle of attack since the errors in lift coefficient 
prediction are directly related to the asymmetry in the boundary layers and the size of the 
wake. 


7.3 Unsteady RANS Analysis 

Given these observations, a few computational strategies can be considered. The most general 
solution would be direct numerical simulation, but the best assessment to date is that a flow 
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Figure 7-6: Time Accurate and Time Averaged Effects of Vortex Origination Point 

field with this complexity could only be run at Reynolds Numbers in the thousands to tens of 
thousands. The next best choice would be a Large Eddy Simulation, but the best assessment 
to date is that a flow field with this complexity could only be reasonably run at Reynolds 
Numbers in the tens to hundreds of thousands. The final potential solution, which is exercised 
here is Unsteady Reynolds Averaged Navier Stokes (URANS). The following quote is a brief 
explanation of the URANS concept from Reference [12]: 

It was noted that RANS models treat all of the unsteadiness in the flow as ‘turbu¬ 
lence’. The concept of unsteady RANS modeling has been proposed; it is sometimes 
claimed that this type of calculation is a sort of very large eddy simulation (VLES). 

Let us discuss this seeming contradiction. Sometimes a RANS calculation fails to 
converge to a steady state solution; in its stead, an unsteady (usually time-periodic) 
solution is produced. Some authors claim that the results can be regarded as a 
VLES but this is hard to justify. If a flow has a separation between a large-scale 
unsteady field and small scale ‘turbulence’ and thus a spectrum with two peaks and a 
gap between them, VLES might make sense. Unfortunately, this kind of flow occurs 
very rarely, if at all. The atmosphere and ocean were thought - behave in this way 
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Figure 7-7: Steady and Unsteady RANS Analysis for B1 Cupped Foil 

but recent evidence suggests that this is not the case. Without a spectral gap it is 
difficult to justify the idea that unsteady RANS calculations are a kind of VLES. 

Also, a model designed to represent only part of the unsteady motion should differ 
from one that models all of the unsteady motion. 

CFDSHIP has a built in time accurate mode of operation that allows URANS calculations. 
The question arises whether or not the flow environment around a foil experiences the type of 
flow mentioned here. The lessons learned in the RANS grid convergence study were used to 
make a grid with higher density in the wake and adequate density on the foil. This grid is 
analyzed using the URANS formulation of CFDSHIP. The vortex shedding patterns observed 
in the calculations show a high degree of asymmetry and a frequency similar to that observed 
in the experiment. A time average of the fluid velocities around the foil is taken and the 
agreement with the experiment is remarkable. Figure 7-7 shows the results of URANS, RANS 
and the experiment. 
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Figure 7-8: Unsteady Rans Non-Isotropic Turbulence 

It should be noted here that the unsteady RANS formulation is very unstable. The actual 
run that is reported took thirty hours to run, but finding a set of run time parameters that 
would allow it to run took months. The horizontal pressure waves reported in Reference [23] 
were observed and seem to be the source of the unstable behavior. 

7.4 Conclusions 

This research indicates that a purely time averaged analysis of the flow field associated with 
a two dimensional foil fails to capture physics which are relevant to determining time aver¬ 
age forces. The experimentally observed and computationally reproduced vortex shedding is 
demonstrated to result in an over-prediction of lift on foils with relevant geometries. 

A series of new experiments and analysis of old experiments yields a correction which sig¬ 
nificantly improves the time average calculations, although it does not fully model the relevant 
physics. The improvement is shown to enhance not only the global lift coefficient predictions, 
but also the local velocity predictions. The experimental time average velocity measurements 
are tabulated to allow future comparisons with flow code improvements. 
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Chapter 8 


Recommendat ions 


8.1 Long Term Solution 

The computational effort required to determine the effects of vortex shedding on a two dimen¬ 
sional foil is considerable. In this research it has been found that approximately 10,000 time 
steps per chord length are required to fully resolve the time accurate vortex shedding patterns. 
For most practical engineering solutions this option is clearly not practical. 

This research indicates that it will be necessary to model anisotropic turbulence in order 
to fully capture the effects of vortex shedding on a two dimensional foil. A steady turbulence 
modeling technique such as second moment closure may be capable of predicting the time 
average flow field correctly since it explicitly solves a set of differential equations for each of the 
turbulence terms. Flow codes that use second moment closure for turbulence modeling are not 
widely available and currently are used primarily in research. 

8.2 Short Term Solution 

Many modern applications use integral boundary layer theory and panel methods to determine 
forces on surfaces that act as a two dimensional foil. For these types of applications Equation 
3.1 can be used to estimate the difference between the steady state viscous calculation and the 
experiment or average time accurate calculation. The following example illustrates the manner 
in which this research can be used: 
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Example (Data from B1 Cupped Foil) 

• A foil is evaluated using a fully viscous IBLT code (such as MSES) to have a lift coefficient 
of 0.64. 

• IBLT gives an upper side displacement thickness of 0.021 chord lengths. 

• IBLT gives a lower side displacement thickness of 0.005 chord lengths. 

1. If the application lends itself to adjusting the computational angle of attack then the 
computational angle of attack is set 0.92 degrees less than the experimental angle of 
attack. This approach was used for all of the comparisons with water tunnel experiments 
reported in this thesis. 


*vortices — 


Kp-SL 

chord 


* 1 
bottom J 


(.021 - .005) 
L0 


= —.016 radians = —0.92 deg 


2. If the application does not lend itself easily to adjusting the foil angle of attack then 
thin foil theory can be used to estimate the change in lift associated with the angle 
error predicted by Equation 3.1. Equation 8.1 shows the approximate change in normal 
coefficient for this case. 


ACx vortices = ~2tt = ~2x(.016) = -0.1 

27r = lift slope from thin foil theory 


Using this method, the lift coefficient (normal vector) is adjusted from 0.64 to 0.54. 
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Appendix A 


Coupled Integral Boundary Layer 
Theory Convergence Study 
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Figure A-l: NACA 66021 Foil Geometry 

A.l Foil Geometry 

The foil used to demonstrate convergence is the NACA 66021 foil. The foil geometry is 
generated based on the mathematical equations in Reference [28]. The individual data points 
experimentally measured are listed in an appendix of the same reference and are shown along 
with the derived surface in Figure A-l. This foil is very thick (21%) and has maximum 
thickness relatively far back on the foil. These characteristics are similar to the shape of 
a propeller blade near the hub and have shown a tendency to result in larger discrepancies 
between experiments and predictions. 


A.2 Grid Convergence 

The domain used during the grid convergence study is eight chord lengths long and 3.75 chord 
lengths high. A base-line coarse mesh is created and then subdivided in streamwise and 
cross stream directions. Figure A-2 shows lift coefficient convergence versus total number 
of nodes. Table A.l shows the final grid dimensions chosen. In reality, the grid chosen is 
far beyond what is required, but processor speeds are such that the excess gridding is of no 
consequence. Processing time for grid generation and analysis ranged from about two seconds 
to 60 seconds. The majority of the calculations require between 5 and 15 Newton iterations to 
reach convergence. 
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Description 

Quantity 

Foil Top and Bottom Nodes 

150 

Upstream Nodes 

50 

Downstream Nodes 

50 

Upper Half of Domain Streamlines 

30 

Lower Half of Domain Streamlines 

30 

Total Number of Nodes 

15000 


Table A.l: Converged Grid Characteristics 
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Figure A-3: Predicted Lift Coefficient vs MSIS Domain Length 

A.3 Domain Size 

The limits of the domain for this study are evaluated by changing the size of the domain and 
observing the effects on the predicted lift. Grid convergence is assured at each new domain 
length. Figures A-3 and A-4 show the variation in predicted lift coefficient based on the domain 
size. A domain height of 3.75 chord lengths is chosen based on the actual tunnel size and a 
distance upstream and downstream of 3.0 chord lengths is chosen based on Figure A-3. 

A.4 Mach Number 

Ultimately the goal of this research is in evaluating incompressible flows relevant to marine 
propellers. In general, incompressible flow is very similar to low Mach Number conditions in 
air, so the same tools and experiments will be useful. The experimental Mach Number for 
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Figure A-4: Predicted Lift Coefficient vs MSIS Domain Height 

the airfoils considered here is typically less than 0.1. Figure A-5 shows that lift coefficient 
predictions are insensitive to Mach Number over the region of interest. The Mach number 
chosen is 0.025 based on anticipated use with incompressible flow. 
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Appendix B 


HRA Fine Trailing Edge Foil 
Experimental Data 
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HRA Fine Trailing Edge Foil Geometry at AOA=0.0 


79 












































































































































































Chord = 405mm 


U/Uinf 


Uinf= 


V/Uinf 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


0.210 


0.111 


- 0.185 


- 0.284 


- 0.383 


- 0.481 


- 0.580 


- 0.679 


- 0.778 


- 0.877 


- 0.877 


- 0.877 


- 0.877 


- 0.877 


- 0.117 


- 0.117 


- 0.117 


- 0.117 


- 0.117 


- 0.110 


- 0.103 


- 0.096 


- 0.090 


- 0.083 


- 0.0761 


- 0.069 


- 0.049 


- 0.042 


- 0.036 


- 0.036 


- 0.0361 


- 0.036 


- 0.036 


- 0.036 


- 0.002 


0.032 


0.066 


0.100 



0.938 


0.921 


0.8701 


0.730 


0.889 


0.958 


0.981 


0.9901 


0.992 


0.994 


0.998 


1.000 


1.0031 


1.005 


- 0.041 


- 0.049 


- 0.060 


0.062 


0.072 


0.075 


0.073 


- 0.066 


0.055 


0.043 


0.0581 


0.026 


0.094 


0.006 


0.162 


0.111 


0.077 


0.056 


0.041 


0.030 


0.0301 


0.032 


0.028 


0.030 


5.312 m/sec 


Cp 


- 0.004 


- 0.006 


- 0.001 


- 0.002 


0.009 


0.024 


0.045 


0.060 


0.068 


0.075 


0.059 


0.076 


0.117 


0.233 


0.092 


0.035 


0.016 


0.008 


0.007 


0.006 


0.002 


0.000 


- 0.004 


- 0.006 


AOA= 


4.3 degrees 


Y/C 


U/Uinf 



- 0.8771 


- 0.778 


- 0.6791 


- 0.580 


- 0.481 


- 0.3831 


- 0.284 


- 0.185 


- 0.0861 


0.012 


0.111 


0.2101 


0.309 


0.407 


0.506 


0.605 


0.704 


0.802 


0.901 


1.0001 


1.099 


1.198 


0.0001 


0.000 


0.100 


0.100 


0.100 


0.100 


0.100 


0.100 


0.094 


0.087 


0.080 


0.073 


0.0661 


0.060 


0.053 


0.046 


0.039 


0.0331 


0.026 


0.0191 


0.019 


0.019 


0.019 


0.019 


0.000 


0.000 



1.252 


1.242 


1.216 


1.186 


1.131 


1.067 


1.0181 


1.009 


1.011 


1.0091 


1.007 


1.005 


0.000 


0.000 


0.030 


0.041 


0.060 


0.0851 


0.126 


0.194 


0.228 


0.098 


0.040 


- 0.049 


- 0.064 


- 0.058 


- 0.1241 


- 0.154 


- 0.181 


- 0.173 


- 0.139 


- 0.094 


- 0.068 


- 0.0561 


- 0.045 


- 0.036 


0.000 


0.000 


Cp 


- 0.006 


- 0.010 


- 0.015 


- 0.016 


- 0.032 


- 0.084 


- 0.268 


- 0.334 


- 0.320 


- 0.318 


- 0.286 


- 0.274 


- 0.247 


- 0.215 


- 0.156 


- 0.085 


- 0.028 


- 0.014 


- 0.013 


- 0.011 


- 0.008 


- 0.006 


0.000 


0.000 


HRA Fine Foil Experimental Data at AOA=4.3 deg 
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Chord = 405mm 


XJC 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.877 


-0.877 


-0.778 


-0.679 


Y/C 


-0.090 


-0.090 


-0.090 


-0.090 


-0.090 


-0.086 


-0.083 


-0.079 


-0.076 


-0.072 


-0.062 


-0.052 


-0.048 


-0.048 


-0.048 


-0.048 


-0.048 


-0.048 


-0.014 


0.020 


0.053 


0.087 


0.087 


0.087 


U/Uinfl 


1.002 


1.004 


1.008 


1.002 


0.993 


0.978 


0.972 


0.968 


0.968 


0.972 


1.017 


0.906 


0.936 


0.974 


0.985 


0.991 


0.995 


0.998 


0.998 


0.998 


1.000 


1.004 


1.002 


1.000 


Uinf=| 


V/Uinf 


-0.025 


-0.030 


-0.036 


-0.038 


-0.040 


-0.032 


-0.023 


- 0.011 


- 0.002 


0.011 


-0.078 


-0.089 


0.038 


0.030 


0.023 


0.017 


0.013 


0.008 


0.011 


0.011 


0.013 


0.015 


0.021 


0.026 


5.288 m/sec 


Cp 


-0.003 


-0.005 


-0.009 


-0.003 


0.006 


0.022 


0.027 


0.031 


0.031 


0.028 


- 0.021 


0.086 


0.061 


0.025 


0.014 


0.009 


0.005 


0.001 


0.001 


0.001 


0.000 


-0.004 


- 0.002 


- 0.001 


AOA= 

1.1 degrees 1 

X/C 

Y/C 

U/Uinf| 


V/Uinf 


-0.679 

0.087 

1.000 

0.026 

-0.580 

0.087 

1.000 

0.038 

-0.481 

0.087 

0.997 

0.059 

-0.383 

0.087 

1.002 

0.098 

-0.284 

0.084 

1.065 

! 6.151 

-0.185 

0.081 

1.137 

0.117 

-0.086 

0.077 

1.163 

j 0.085 

0.012 

0.074 

1.184 

0.017 

0.111 

0.070 

1.178 

0.002 

0.210 

0.067 

1.180 

-0.009 

0.309 

0.063 

1.171 

-0.061 

0.407 

0.060 

1.157 

-0.096 

0.506 

0.056 

1.114 

-0.132 

0.605 

0.053 

1.053 

-0.136 

0.704 

0.050 

1.012 

-0.100 



-0.061 


-0.042 


-0.032 


-0.026 


-0.019 


Cp 


- 0.001 


- 0.001 


0.002 


-0.007 


-0.078 


-0.153 


-0.180 


- 0.201 


-0.194 


-0.196 


-0.187 


-0.174 


-0.129 


-0.064 


-0.017 


-0.008 


-0.009 


-0.008 


-0.006 


- 0.001 


HRA Fine Foil Experimental Data at AOA=l.l deg 
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x/c 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


0.309 


0.111 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.877 


-0.877 


Y/CI 


-0.034 


-0.034 


-0.034 


-0.0341 


-0.034 


-0.037 


-0.041 


-0.044 


-0.0481 


-0.051 


-0.051 


-0.058 


-0.068 


-0.0721 


-0.075 


-0.075 


-0.075 


-0.075 


-0.0751 


-0.075 


-0.041 


-0.0071 


0.027 


0.061 


U/Uin 


0.9861 


0.994 


0.999 


1.001 


0.992 


0.9801 


0.984 


0.992 


1.0031 


1.005 


1.005 


1.048 


1.160 


1.086 


0.984 


0.990 


0.995 


0.9991 


1.001 


0.9991 


1.001 


1.001 


0.9991 


0.999 


V/Uinf| 


-0.028 


-0.026 


-0.0301 


-0.028 


-0.025 


-0.013 


0.011 


0.025 


0.0341 


0.051 


0.051 


0.074 


-0.0571 


-0.176 


- 0.100 


-0.049 


-0.034 


-0.0361 


-0.030 


-0.023 


-0.025 


-0.017 


-0.0151 


- 0.011 


Cp 


0.014 


0.006 


0.000 


- 0.002 


0.008 


0.019 


0.016 


0.008 


-0.004 


-0.006 


-0.006 


-0.052 


-0.175 


-0.106 


0.011 


0.009 


0.004 


0.000 


- 0.002 


0.000 


- 0.001 


- 0.001 


0.001 


0.001 


xq 


-0.778 


-0.6791 


-0.580 


-0.481 


-0.383 


-0.284 


-0.1851 


0.012 


0.111 


0.210 


0.309 


0.4071 


0.506 


0.605 


0.7041 


0.802 


0.901 


1.0001 


1.099 


1.1981 


0.061 


0.061 


0.061 


0.061 


0.061 


0.0641 


0.068 


0.074 


0.078 


0.081 


0.085 


0.0881 


0.092 


0.095 


0.099 


0.1021 


0.102 


0.102 


0.102 


0.1021 


U/Uinf| 


0.999 


0.999 


0.990 


0.980 


0.9461 


0.939 


1.014 


1.100 


1.107 


1.1171 


1.118 


1.118 


1.088 


1.0431 


1.009 


1.0071 


1.009 


1.0091 


1.007 


1.001 


V/Uin 


- 0.011 


- 0.011 


-0.0171 


-0.013 


- 0.002 


0.097 


0.1271 


0.074 


0.053 


0.034 


-0.0041 


-0.044 


-0.079 


-0.085 


-0.0571 


-0.032 


-0.015 


- 0.011 


-0.006 


-0.0081 


Cp 


0.001 


0.001 


0.010 


0.019 


0.052 


0.055 


-0.023 


-0.107 


-0.114 


-0.124 


-0.125 


-0.126 


-0.095 


-0.047 


- 0.010 


-0.007 


-0.009 


-0.009 


-0.007 


- 0.001 



HRA Fine Foil Experimental Data at AOA=-1.7 deg 
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Appendix C 


HRA Blunt Trailing Edge Foil 
Experimental Data 
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Y/C=0.625 



HRA Blunt Foil Geometry at AOA=0.0 
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0.210 


0.111 


-0.086 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.094 


-0.081 


-0.075 


-0.068 


-0.062 


-0.062 


-0.062 


-0.062 


-0.062 


-0.062 


-0.019 


0.025 


U/Uinfl 


0.998 




















0.938 


0.932 


0.924 


0.873 


0.789 


0.879 


0.943 


0.9691 


0.981 


0.990 


0.994 


0.994| 


1.002 


-0.045 


-0.039 


-0.078 


-0.078 


0.006 


0.117 


0.097 


0.074 


0.055 


0.041 


0.031 


0.035 


0.037 


0.055 


-0.037 


-0.051 


-0.090 


-0.115 


-0.142 


-0.172 


-0.164 


- 0.111 


-0.074 


-0.055 


-0.045 


-0.037 


Cp 


-0.003 


-0.005 


-0.007 


- 0.011 


-0.014 


- 0.020 


-0.033 


-0.080 


-0.231 


-0.306 



-0.299 


-0.264 


-0.253 


-0.232 


- 0.211 


-0.185 


-0.134 


-0.059 


-0.014 


-0.013 


-0.015 


-0.013 


- 0.012 


HRA Blunt Foil Experimental Data at AOA=3.6 deg 































































































































































































































































HRA Blunt Foil Experimental Data at AOA=0.9 
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Chord = | 


x/c 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


0.210 


0.111 


-0.086 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


405mm 


Y/C 


-0.074 


-0.074 


-0.074 


-0.074 


-0.074 


-0.076 


-0.077 


-0.079 


-0.080 


-0.082 


-0.083 


-0.085 


-0.088 


-0.089 


-0.091 


-0.092 


-0.092 


-0.092 


Uinf=|5. 


V/Uinf 


13 m/sec 


Cp 


AOA= 


y/ c 


-0.877 


-0.877 


-0.778 


-0.679 


-0.580 


-0.481 


-0.383 


-0.284 


-0.185 


-0.086 


0.012 


0.111 


0.210 


0.309 


0.407 


0.506 


0.605 


0.704 


0.802 


0.901 


1.000 


1.099 


0.9 degrees 


Y/C U/Uinf 


V/Uinf 

Cp 

-0.004 

0.000 

-0.002 

-0.002 

0.000 

0.002 

0.000 

0.002 

0.004 

0.011 

0.016 

0.019 

0.045 

0.033 

0.113 

0.007 

0.123 

-0.062 

0.105 

-0.098 

0.045 

-0.119 

0.037 

-0.127 

0.023 

-0.136 

-0.010 

-0.140 

-0.037 

-0.140 

-0.076 

-0.134 

-0.117 

-0.105 

-0.125 

-0.044 

-0.080 

-0.006 

-0.049 

-0.005 

-0.035 

-0.011 

-0.029 

-0.013 


U/Uinf 


1.000 

-0.021 

-0.001 

1.004 

-0.023 

-0.005 

0.995 

-0.020 

0.005 

0.991 

-0.018 

0.009 

0.981 

0.000 

0.019 

0.995 

0.031 

0.005 

1.004 

0.018 

-0.004 

1.008 

0.016 

-0.008 

1.014 

0.021 

-0.014 

1.020 

0.029 

-0.021 

1.037 

0.008 

-0.038 

1.045 

0.033 

-0.047 

1.088 

-0.029 

-0.092 

1.086 

-0.080 

-0.093 

1.014 

-0.117 

-0.021 

0.971 

-0.053 

0.027 

0.981 

-0.023 

0.019 

0.991 

-0.012 

0.009 

0.995 

-0.012 

0.005 

0.998 

-0.012 

0.001 

0.998 

-0.010 

0.002 

0.998 

-0.002 

0.002 

1.000 

-0.004 

0.000 


0.037| 


0.081 


0.081 


0.081 


0.081 


0.081 


0.081 


0.082 


0.084 


0.085 


0.087 


0.088 


0.090 


0.091 


0.093 


0.094 


0.095 


0.097 


0.098 


0.098 


0.098| 


0.098 


1.000 


1.002 


0.998 


0.998 


0.989 


0.981 


0.965 


0.987 


1.053 


1.088 


1.112 


1.119 


1.127 


1.131 


1.131 


1.123 


1.094 


1.035 


1.002 


1.004 


1.010 


1.012 


HRA Blunt Foil Experimental Data at AOA=-0.9 
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Appendix D 


HRA Beveled Trailing Edge Foil 
Experimental Data 
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0.7000 


0.6950 


0.6899 


0.6849 


0.6748 


0.6748 


0.6648 


0.6498 


0.6248 


0.5997 


0.5747 


0.5497 


0.5247 


0.4997 


0.4747 


0.4496 


0.4246 


0.3996 


0.3496 


0.2996 


0.2496 


0.1996 


0.1496 


0.0996 


0.0496 


-0.0003 


-0.0503 


-0.1003 


-0.0075 


-0.00331 


0.0009 


0.0051 


0.0136 


0.0136 


0.0157 


0.0187 


0.0231 


0.0267 


0.0298 


0.0323| 


0.0344 


0.0363 


0.0379 


0.0392 


0.0404 


0.0414 


0.0428 


0.0436 


0.0438 


0.0435 


0.0427 


0.0414 


0.0394| 


0.0368 


0.0336 


0.0294 


-0.1253 


-0.1503 


-0.1752 


- 0.2002 


-0.2252 


-0.2501 


-0.2651 


-0.2751 


-0.2851 


-0.2901 


-0.2950 


-0.2975 


-0.2990 


-0.3000 


-0.2990 


-0.2975 


-0.2950 


-0.2899 


-0.2849 


-0.2749 


-0.2649 


-0.2499 


-0.2248 


-0.1998 


-0.1748 


-0.1498 


-0.1248 


-0.0998 


0.0270 


0.0243 


0.0212 


0.0178 


0.0138 


0.0090 


0.0056 


0.0029 


-0.0003 


- 0.0022 


-0.0046 


-0.0063 


-0.0076 


-0.0099 


- 0.0121 


-0.0134 


-0.0149 


-0.0171 


-0.0188 


-0.0215 


-0.0237 


-0.0264 


-0.0299 


-0.0325 


-0.0343 


-0.0356 


-0.0363 


-0.0365 


-0.0498 


0.0002 


0.0502 


0.1002 


0.1502 


0.2001 


0.2501 


0.3001 


0.3501 


0.4001 


0.4251 


0.4501 


0.4751 


0.5001 


0.5251 


0.5501 


0.5751 


0.6001 


0.6251 


0.6501 


0.6651 


0.6751 


0.6851 


0.6900 


0.6950 


0.7000 


Y/C 


-0.0359 


-0.0345 


-0.0325 


-0.0303 


-0.0280 


-0.0257 


-0.0236 


-0.0217 


- 0.0200 


-0.0186 


-0.0180 


-0.0175 


-0.0169 


-0.0164 


-0.0160 


-0.0155 


-0.0151 


-0.0146 


-0.0141 


-0.0135 


-0.0131 


-0.0128 


-0.0126 


-0.0124 


-0.0123 


- 0.0122 
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Chord = 1405mm 


X/CI Y/C 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


0.210 


0.111 


0.074 


-0.086 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.142 


-0.142 


-0.142 


-0.142 


-0.142 


-0.135 


-0.128 


- 0.122 


-0.115 


-0.108 


- 0.101 


-0.094 


-0.092 


-0.080 


-0.074 


-0.067 


-0.060 


-0.060 


-0.060 


-0.060 


-0.060 


-0.060 


-0.017 


0.027 


U/Uin 


0.997 


1.001 


0.9931 


0.985 


0.974 


0.984 


0.978 


0.964 


0.9481 


0.940 


0.938 


0.934 


0.934 


0.925 


0.874 


0.783 


0.878 


0.944 


0.972 


0.982 


0.9871 


0.991 


0.9931 


0.999 


Uinf= 5. 


V/Uinf 


094 m/sec 


-0.0391 


-0.047 


-0.049 


-0.049 


-0.043 


-0.0351 


-0.053 


-0.057 


-0.057 


-0.049 


-0.047 


-0.041 


-0.0331 


-0.079 


-0.082 


0.004 


0.1221 


0.098 


0.073 


0.0571 


0.043 


0.033 


0.035 


0.035 


0.002 


- 0.002 


0.005 


0.013 


0.025 


0.016 


0.021 


0.034 


0.049 


0.057 


0.059 


0.063 


0.063 


0.069 


0.115 


0.193 


0.108 


0.049 


0.025 


0.017 


0.012 


0.008 


0.006 


0.000 


I 


-0.8771 


-0.877 


-0.877 


-0.778 


-0.679 


-0.580 


-0.481 


-0.3831 


-0.284 


-0.185 


-0.086 


0.074 


0.111 


0.2101 


0.309 


0.407 


0.481 


0.6301 


0.704 


0.802 


0.901 


1.0001 


1.099 


1.198 


0.027 


0.0701 


0.113 


0.113 


0.1131 


0.113 


0.113 


0.113 


0.106 


0.099 


0.092 


0.081 


0.079 


0.072 


0.065 


0.058 


0.053 


0.043 


0.038 


0.031 


0.031 


0.031 


0.031 


0.031 


U/Uinfl 


0.999 


1.005 


1.011 


1.011 


1.013 


1.0171 


1.029 


1.066 


1.197 


1.270 


V/Uinf| 


0.035 


0.033 


0.035 


0.0471 


0.065 


0.088 


0.1241 


0.183 


0.206 


0.100 


0.043 


-0.049 


-0.051 


-0.0631 


- 0.102 


-0.124 


-0.145 


-0.179 


-0.165 


-0.116 


-0.077 


-0.057 


-0.045 


-0.037 


Cp 


0.000 


-0.006 


- 0.012 


- 0.012 


-0.015 


- 0.021 


-0.037 


-0.085 


-0.238 


-0.312 


-0.305 


-0.278 


-0.269 


-0.255 


-0.236 


-0.213 


-0.193 


- 0.112 


-0.059 


-0.016 


- 0.012 


-0.015 


-0.016 


- 0.010 


HRA Beveled Foil Experimental Data at AOA=4.4 deg 
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HRA Beveled Foil Experimental Data at AOA=3.35 deg 












































































































































































































































































Chord = 


405mm 


Y/C 


1.198 


1.099 


1.000 


0.901 


0.802 


0.704 


0.605 


0.506 


0.407 


0.309 


0.210 


0.111 


0.074 


-0.086 


-0.185 


-0.284 


-0.383 


-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.085 


-0.085 


-0.085 


-0.085 


-0.0851 


-0.085 


-0.086 


-0.086 


-0.086 


-0.086 


-0.086 


-0.086 


-0.086 


-0.087 


-0.087 


-0.087 


-0.087 


-0.0871 


-0.087 


-0.087 


-0.087 


-0.087 


-0.044 


- 0.001 


U/Uinf| 


0.994 


1.000 


0.996 


0.988 


0.976 


0.990 


1.002 


1.0041 


1.004 


1.008 


1.024 


1.031 


1.037 


1.0671 


1.057 


0.986 


0.955 


0.976 


0.988 


0.994 


0.996 


0.998 


0.9981 


1.000 


Uinf= 


V/Uinf 


-0.026 


-0.024 


-0.026 


-0.024 


-0.006 


0.016 


0.010 


0.010 


0.014 


0.022 


0.006 


0.026 


0.033 


-0.0261 


-0.063 


-0.091 


- 0.022 


- 0.002 


0.006 


0.006 


- 0.002 


- 0.002 


0.0001 


0.002 


5.08 m/sec 


Cp 


0.006 


0.000 


0.004 


0.011 


0.023 


0.010 


- 0.002 


-0.004 


-0.004 


-0.008 


-0.024 


-0.032 


-0.039 


-0.069 


-0.061 


0.010 


0.044 


0.023 


0.012 


0.006 


0.004 


0.002 


0.002 


0.000 


AOA= 


yj c 


-0.8771 


-0.877 


-0.877 


-0.778 


-0.679 


-0.5801 


-0.481 


-0.3831 


-0.284 


-0.185 


-0.086 


0.074 


0.111 


0.2101 


0.309 


0.407 


0.481 


0.6301 


0.704 


0.802 


0.901 


1.0001 


1.099 


1.1981 


0.4 degrees 


Y/C 


- 0.001 


0.0431 


0.086 


0.086 


0.086 


0.086 


0.0861 


0.086 


0.086 


0.0861 


0.086 


0.0871 


0.087 


0.087 


0.087 


0.087 


0.087 


0.0871 


0.087 


0.088 


0.088 


0.088 


0.088 


0.088 


1.000 


1.0001 


1.004 


1.002 


1.000 


0.996 


0.988 


0.980 


1.018 


1.081 


1.112 


1.136 


1.138 


1.1421 


1.144 


1.142 


1.136 


1.087 


1.041 


1.008 


1.0121 


1.014 


1.014 


1.0101 


0.002 


0.004 


0.0081 


0.010 


0.014 


0.022 


0.035 


0.063 


0.124 


0.116 


0.094 


0.0301 


0.028 


0.008 


-0.018 


-0.043 


-0.067 


- 0.120 


-0.116 


-0.071 


-0.041 


-0.0261 


- 0.020 


-0.014 


Cp 


0.000 


0.000 


-0.004 


- 0.002 


0.000 


0.004 


0.011 


0.018 


-0.026 


-0.091 


-0.123 


-0.145 


-0.148 


-0.152 


-0.154 


-0.153 


-0.147 


-0.098 


-0.049 


- 0.010 


-0.013 


•0.014 


-0.014 


-0.01Q 


HRA Beveled Foil Experimental Data at AOA=0.4 deg 
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Chord = I 


y/c 


405mm 


Y/C 



-0.481 


-0.580 


-0.679 


-0.778 


-0.877 


-0.877 


-0.877 


-0.877 


HRA Beveled Foil Experimental Data at AOA=-1.9 deg 
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Appendix E 


B1 Blunt Trailing Edge Foil 
Experimental Data 
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0.6499 


0.6361 


0.6224 


0.6086 


0.5948 


0.5810 


0.5672 


0.5534 


0.5374 


0.5060 


0.4745 


0.4436 


0.4140 


0.3853 


0.3571 


0.3292 


0.3010 


0.2721 


0.2421 


0.2107 


0.1778 


0.1437 


0.1088 


0.0731 


0.0372 


0.0011 


-0.0349 


-0.0703 


Y/C 


- 0.0212 


-0.0228 


-0.0241 


-0.0253 


-0.0263 


-0.0273 


-0.0281 


-0.0288 


-0.0296 


-0.0309 


-0.0320 


-0.0329 


-0.0338 


-0.0346 


-0.0353 


-0.0360 


-0.0365 


-0.0370 


-0.0374 


-0.0378 


-0.0380 


-0.0382 


-0.0382 


-0.0382 


-0.0381 


-0.0378 


-0.0374 


-0.0369 


X/C 


-0.1046 


-0.1374 


-0.1681 


-0.1963 


-0.2215 


-0.2432 


-0.2612 


-0.2762 


-0.2884 


-0.2984 


-0.3066 


-0.3136 


-0.3197 


-0.3252 


-0.3303 


-0.3348 


-0.3387 


-0.3422 


-0.3450 


-0.3473 


-0.3490 


-0.3499 


-0.3499 


-0.3489 


-0.3468 


-0.3433 


-0.3383 


0.3317 


Y/C 


-0.0364 


-0.0357 


-0.0350 


-0.0343 


-0.0336 


-0.0328 


-0.0321 


-0.0314 


-0.0307 


-0.0301 


-0.0294 


-0.0287 


-0.0281 


-0.0275 


-0.0268 


-0.0260 


-0.0252 


-0.0242 


-0.0231 


-0.0218 


-0.0203 


-0.0185 


-0.0165 


-0.0144 


- 0.0120 


-0.0093 


-0.0065 


0.0034 


X/C 


-0.3233 


-0.3130 


-0.3006 


-0.2859 


-0.2688 


-0.2490 


-0.2266 


-0.2016 


-0.1744 


-0.1453 


-0.1147 


-0.0829 


-0.0502 


-0.0168 


0.0169 


0.0507 


0.0844 


0.1179 


0.1509 


0.1832 


0.2147 


0.2451 


0.2745 


0.3030 


0.3304 


0.3567 


0.3820 


0.4062 


- 0.0002 


0.0032 


0.0068| 


0.0105 


0.01431 


0.0182 


0.0221 


0.0260 


0.0298 


0.0335| 


0.0369 


0.0400 


0.0428 


0.0451 


0.04701 


0.0484 


0.0495 


0.0501 


0.0503 


0.0501 


0.0496 


0.0487 


0.0475 


0.0459 


0.0440 


0.0418 


0.0392 


0.0363 


0.4293 


0.4511 


0.4715 


0.4905 


0.5080 


0.5237 


0.5376 


0.5499 


0.5568 


0.5706 


0.5844 


0.5982 


0.6120 


0.6258 


0.6396 


0.6499 


Y/C 


0.0332 


0.0299 


0.0265 


0.0231 


0.0198 


0.0166 


0.0137 


0.0111 


0.0095 


0.0063 


0.0030 


-0.0004 


-0.0040 


-0.0076 


- 0.0112 


-0.0140 



B1 Blunt Foil Geometry at AOA=0.0 

































































































































































































































Chord = 


y/c 


455 mm 


Y/C 


Uinf= 5. 


V/Uin 


114 m/sec 


AOA= 


yjc 


5.2 degrees 


U/Uinf 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.780 


-0.136 


-0.136 


-0.1361 


0.136 


0.136 


0.131 


0.125 


0.119 


0.1131 


0.107 


0.101 


0.095 


0.083 


0.077 


0.071 


0.0651 


0.065 


0.065 


0.065 


0.065 


0.065 


0.021 


0.023 


0.067 


0.9971 


1.005 


0.980 


0.980 


0.962 


0.980 


0.993 


0.987 


0.974 


0.960 


0.043 


0.045 


0.055 


0.049 


0.037 


0.027 


0.049 


0.059 


0.0631 


0.070 


0.072 


0.078 


0.096 


0.078 


0.057 


0.041 


0.1371 


0.117 


0.092 


0.065 


0.055 


0.061 


0.065 


0.061 


0.002 


-0.006 


0.019 


0.019 


0.037 


0.020 


0.005 


0.011 


0.024 


0.037 


0.048 


0.062 


0.101 


0.125 


0.163 


0.207 


0.127 


0.060 


0.033 


0.020 


0.015 


0.007 


- 0.001 


-0.013 


-0.780 


-0.692 


-0.6041 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 



0.111 


0.111 


0.111 


0.111 


0.111 


0.111 


0.1051 


0.099 


0.093 


0.087 


0.075 


0.069 


0.063 


0.057 


0.051 


0.045 


0.039 


0.039 


0.039 


0.039 


0.039 


U/Uinf 


1.015 


1.021 


1.0271 


1.038 


1.089 


1.240 


1.332 


1.345 


1.353 


1.361 


1.302 


1.263 


1.218 


1.156 


1.087 


1.032 


1.013 


1.007 


1.009 


1.009 


0.997 


V/Uinf| 


0.065 


0.084 


0.113 


0.162 


0.2331 


0.258 


0.158 


0.078 


0.029 


-0.051 


-0.108 


-0.160 


-0.196 


-0.215 


-0.209 


-0.1741 


-0.131 


-0.092 


-0.0701 


-0.061 


-0.0491 


Op 


-0.017 


-0.024 


-0.033 


-0.052 


- 0.120 


-0.302 


-0.399 


-0.408 


-0.416 


-0.427 


-0.354 


-0.311 


-0.261 


-0.191 


-0.113 


-0.048 


- 0.022 


- 0.011 


- 0.012 


- 0.011 


0.002 


B1 Blunt Foil Experimental Data at AOA=5.2 deg 
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Chord =j 


X/C 


455 mm 


Y/C 


Uinf= 5. 


V/Uinf 


26 m/sec 


.25 degrees 


Y/C U/Uinf 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.108 


-0.108 


-0.108 


-0.108 


-0.108 


-0.105 


-0.103 


- 0.101 


-0.098 


-0.096 


-0.093 


-0.091 


-0.088 


-0.086 


-0.0831 


-0.081 


-0.078 


-0.078 


-0.078 


-0.078 


-0.078 


-0.078 


-0.034 


0.010 


U/Uinf! 


0.998 


1.004 


0.998| 


0.989 


0.975 


0.987 


1.004 


1.010 


1.004 


1.000 


1.002 


0.992 


0.994 


0.979 


0.960 


0.935 


0.892 


0.914 


0.958| 


0.977 


0.989 


0.992 


0.996 


0.998 


-0.051 


-0.053 


-0.053 


-0.048 


-0.032 


-0.008 


-0.017 


-0.023 


-0.032 


-0.040 


-0.046| 


-0.051 


-0.048 


-0.074 


-0.070 


-0.070 


-0.046 


0.034 


0.044 


0.036| 


0.029 


0.017 


0.025 


0.034 


0.001 


-0.005 


0.000 


0.010 


0.024 


0.013 


-0.004 


- 0.010 


-0.004 


- 0.001 


-0.003 


0.006 


0.005 


0.018 


0.037 


0.060 


0.101 


0.081 


0.040 


0.022 


0.011 


0.007 


0.003 


0.001 


-0.780 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187| 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


0.053 


0.097 


0.097 


0.097 


0.097| 


0.097 


0.097 


0.095 


0.093 


0.090| 


0.088 


0.085 


0.083 


0.080 


0.078 


0.075 


0.073 


0.070 


0.068 


0.068 


0.068 


0.068 


0.068 


1.004 


1.010 


1.010 


1.010 


1.013 


1.029 


1.110 


1.198 


1.230 


1.249 



V/Uinf 


0.032 


0.036 


0.048 


0.065 


0.097 


0.148 


0.196 


0.158 


0.105 


0.068 


0.002 


- 0.021 


-0.061 


-0.105 


-0.148 


-0.1771 


-0.181 


-0.148 


-0.106 


-0.074 


-0.0531 


-0.046 


-0.040 


Cp 


-0.004 


- 0.010 


- 0.011 


- 0.012 


-0.018 


-0.040 


-0.136 


-0.230 


-0.262 


-0.282 


-0.294 


-0.304 


-0.287 


-0.274 


-0.244 


-0.188 


- 0.112 


-0.050 


-0.019 


-0.014 


- 0.021 


- 0.020 


-0.014 


B1 Blunt Foil Experimental Data at AOA=2.25 deg 
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Chord = 


x/c 


455 mm 


Y/C 


U/Uinfl 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


- 0.101 


- 0.101 


- 0.101 


- 0.101 


- 0.101 


-0.099 


-0.098 


-0.096 


-0.094 


-0.093 


-0.091 


-0.0901 


-0.088 


-0.087 


-0.085 


-0.083 


-0.082 


-0.082 


-0.082 


-0.082 


-0.082 


-0.082 


-0.038 


0.006 



1.039 


1.031 


1.0431 


1.035 


1.029 


1.0241 


0.985 


0.956 


0.973 


0.985 


0.993 


0.997 


0.998 


1.000 


Uinf= 


V/Uin 


-0.041 


-0.037 


-0.035 


-0.027 


- 0.012 


0.012 


0.014 


0.008 


- 0.002 


- 0.002 


-0.014 


-0.015 


-0.017 


-0.042 


-0.042 


-0.056 


-0.079 


-0.027 


0.000 


0.002 


0.002 


0.000 


0.006 


0.012 


5.178 m/sec 


Cp 


-0.007 


-0.009 


-0.003 


0.003 


0.015 


0.003 


-0.018 


-0.028 


-0.032 


-0.036 


-0.040 


-0.032 


-0.044 


-0.037 


-0.031 


-0.025 


0.012 


0.043 


0.026 


0.015 


0.007 


0.003 


0.002 


0.000 


AOA= 0. 


yjc 


4 degrees 


Y 


I 


-0.780 


-0.780 


-0.692 


-0.6041 


-0.516 


-0.429 


-0.341 


-0.2531 


-0.165 


-0.077 


0.011 


0.0991 


0.187 


0.275 


0.363 


0.451 


0.538 


0.6261 


0.714 


0.802 


0.8901 


0.978 


1.066 


0.0501 


0.094 


0.094 


0.094 


0.094 


0.094 


0.094 


0.092 


0.091 


0.0891 


0.088 


0.086 


0.085 


0.083 


0.081 


0.0801 


0.078 


0.077 


0.075 


0.075 


0.075 


0.075 


0.075 


1.000 


1.004 


1.002 


0.998 


0.9931 


0.989 


1.024 


1.099 


1.1471 


1.178 


1.199 


1.2131 


1.211 


1.2071 


1.188 


1.143 


1.080 


1.022 


1.002 


1.0041 


1.010 


1.010 


1.006 


0.012 


0.017 


0.021 


0.031 


0.048 


0.0871 


0.147 


0.153 


0.126 


0.098 


0.041 


0.021 


-0.0151 


-0.060 


- 0.100 


-0.145 


-0.154 


-0.1291 


-0.081 


-0.0501 


-0.037 


-0.029 


-0.025 


Cp 


0.000 


-0.004 


-0.003 


0.001 


0.006 


0.007 


-0.035 


-0.115 


-0.166 


-0.199 


- 0.220 


-0.236 


-0.233 


-0.230 


- 0.210 


-0.164 


-0.095 


-0.030 


-0.006 


-0.006 


- 0.011 


- 0.011 


-0.007 


B1 Blunt Foil Experimental Data at AOA=0.4 deg 
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- 1.5 - 1.25 -1 - 0.75 - 0.5 - 0.25 9 


Q.- 0 - 2 ' 

O rv 



Chord = | 


YJ C 


455 mm 


Y/C 


Uinf=j5. 


V/Uinf 


268 m/sec 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.075 


-0.075 


-0.075 


-0.075 


-0.075 


-0.077 


-0.078 


-0.080 


-0.081 


-0.083 


-0.085 


-0.086 


-0.088 


-0.089 


-0.091 


-0.093 


-0.094 


-0.094 


-0.094 


-0.094 


-0.094 


-0.094 


-0.050 


-0.006 


U/Uinfl 


1.004 


1.006 


1.006 


0.995 


0.987 


1.000 



1.109 


1.002 


0.995 


0.998 


0.998 


1.000 


1.002 


1.000 


-0.023 


- 0.021 


-0.013 


-0.006 


0.013 


0.051 


0.055 


0.049 


0.049 


0.042 


0.0281 


0.030 


0.032 


- 0.002 


-0.006 


-0.038 


-0.129 


-0.103 


-0.057 


-0.032 


-0.023 


-0.017 


-0.015 


-0.008 



- 0.002 


-0.027 


-0.046 


-0.060 


-0.070 


-0.086 


-0.090 


-0.105 


-0.108 


-0.119 


-0.143 


-0.123 


-0.008 


0.004 


0.001 


0.001 


- 0.001 


- 0.002 


0.000 


X/C 


-0.780 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.1651 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


Y/C 


0.038 


0.082 


0.082 


0.082 


0.082 


0.082 


0.082 


0.083 


0.085 


0.087 


0.088 


0.090 


0.091 


0.093 


0.094 


0.096 


0.098 


0.099 


0.101 


0.101 


0.101 


0.101 


0.101 


U/Uinf 


0.998 


0.998 


0.995 


0.989 


0.972 


0.943 


0.940 


1.004 


1.063 


1.109 


V/Uinf 


-0.006 


- 0.002 


- 0.002 


- 0.002 


0.004 



Cp 


0.002 


0.002 


0.005 


0.011 


0.028 


0.055 


0.053 


-0.015 


-0.076 


-0.123 


-0.156 


-0.172 


-0.184 


-0.193 


-0.184 


-0.146 


-0.086 


-0.029 


-0.008 


-0.009 


-0.014 


-0.016 


- 0.010 


B1 Blunt Foil Experimental Data at AOA=-1.95 deg 
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Appendix F 


B1 Antisinging Trailing Edge Foil 
Experimental Data 


100 


0.6500 


0.6435 


0.6367 


0.6297 


0.6226 


0.6153 


0.6080 


0.6007 


0.5935 


0.5861 


0.5784 


0.5699 


0.5605 


0.5499 


0.5377 


0.5238 


0.5080 


0.4906 


0.4716 


0.4511 


0.4293 


0.4062 


0.3820 


0.3568 


0.3304 


0.3030 


0.2746 


-0.0193 


-0.0164 


-0.0135 


-0.0107 


-0.0079 


-0.0053 


-0.0028 


-0.0004 


0.0018 


0.0039 


0.0060 


0.0081 


0.0103 


0.0126 


0.0152 


0.0181 


0.0212 


0.0244 


0.0278 


0.0311 


0.0344 


0.0374 


0.0402 


0.0427 


0.0448 


0.0466 


0.0481 


X/C 


0.2147 


0.1833 


0.1509 


0.1179 


0.0845 


0.0507 


0.0169 


-0.0168 


-0.0501 


-0.0829 


-0.1147 


-0.1453 


-0.1744 


-0.2016 


-0.2265 


-0.2490 


-0.2688 


-0.2859 


-0.3006 


-0.3130 


-0.3233 


-0.3317 


-0.3383 


-0.3433 


-0.3468 


-0.3489 


-0.3499 


Y/C 


0.0500 


0.0505 


0.0505 


0.0502 


0.0495 


0.0484 


0.0468 


0.0448 


0.0424 


0.0395 


0.0363 


0.0328 


0.0290 


0.0251 


0.0211 


0.0171 


0.0132 


0.0094 


0.0056 


0.0020 


-0.0015 


-0.0047 


-0.0078 


-0.0107 


-0.0133 


-0.0157 


-0.0179 


X/C 


-0.3490 


-0.3473 


-0.3450 


-0.3422 


-0.3387 


-0.3348 


-0.3303 


-0.3252 


-0.3197 


-0.3136 


-0.3066 


-0.2984 


-0.2884 


-0.2761 


-0.2612 


-0.2431 


-0.2215 


-0.1963 


-0.1681 


-0.1374 


-0.1046 


-0.0702 


-0.0349 


0.0011 


0.0372 


0.0732 


0.1088 


Y/C 


-0.0216 


-0.0231 


-0.0244 


-0.0256 


-0.0265 


-0.0273 


-0.0281 


-0.0287 


-0.0294 


-0.0300 


-0.0306 


-0.0313 


-0.0319 


-0.0325 


-0.0332 


-0.0338 


-0.0345 


-0.0352 


-0.0358 


-0.0364 


-0.0369 


-0.0374 


-0.0378 


-0.0380 


-0.0382 


-0.0382 


-0.0381 


0.1778 


0.2107 


0.2422 


0.2722 


0.3010 


0.3292 


0.3572 


0.3853 


0.4140 


0.4437 


0.4745 


0.5060 


0.5375 


0.5682 


0.5977 


0.6252 


0.6500 


Y/C 


-0.0377 


-0.0373 


-0.0369 


-0.0364 


-0.0358 


-0.0351 


-0.0344 


-0.0336 


-0.0327 


-0.0317 


-0.0307 


-0.0295 


-0.0281 


-0.0264 


-0.0245 


- 0.0221 


-0.0193 



0.24511 0.0493 -0.3499|| -0.0199| 0.1438| 


B1 Antisinging Foil Geometry at AOA=0.0 





















































































































































































































Chord = 455 mm 


YJ C Y/C 


Uinf= 5 


V/Uinf 


.192 m/sec 


AOA= 5. 


YJC 


25 degrees 


Y/C U/Uin 



1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.1361 


-0.136 


-0.136 


-0.136 


-0.136 


-0.131 


-0.125 


-0.119 


-0.113 


-0.107 


- 0.101 


-0.095 


-0.089 


-0.083 


-0.077 


-0.071 


-0.065 


-0.065 


-0.065 


-0.065 


-0.065 


-0.065 


- 0.021 


0.023 


U/Uinf| 


1.000 


0.998 


0.994| 


0.982 


0.973 


0.992 


1.000 


0.984 


0.973 


0.957 


0.942 


0.924 


0.921 


0.886 


0.855 


0.813| 


0.759 


0.855 


0.928 


0.963| 


0.977 


0.986 


0.990 


1.000 


-0.050 


-0.050 


-0.048 


-0.046 


-0.033 


-0.023 


-0.048 


-0.064 


-0.064 


-0.075 


-0.083 


-0.081 


-0.075 


-0.0961 


-0.079 


-0.060 


0.037 


0.133 


0.112 


0.087 


0.062 


0.046 


0.056 


0.056 


- 0.001 


0.001 


0.005 


0.016 


0.026 


0.008 


- 0.001 


0.014 


0.025 


0.039 


0.053 


0.069 


0.073 


0.103 


0.131 


0.168 


0.211 


0.126 


0.063 


0.033 


0.021 


0.013 


0.008 


- 0.001 


-0.7801 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363| 


0.451 


0.538 


0.626| 


0.714 


0.802 


0.890| 


0.978 


1.066 


0.067 


0.111 


0.111 


0.111 


0.111 


0.111 


0.111 


0.105 


0.099 


0.093| 


0.087 


0.081 


0.075 


0.069 


0.0631 


0.057| 


0.051 


0.045 


0.039 


0.039| 


0.039 


0.039 


0.039 


1.007 


1.017 


1.021 


1.027 


1.042 


1.088 


1.238 


1.329| 


1.341 


1.341 


1.333 


1.333 


1.292 


1.267 


1.225 


1.161 


1.088 


1.040 


1.017 


1.015 


1.017 


1.017 


1.013 


0.060 


0.060 


0.083 


0.112 


0.158 


0.227 


0.248 


0.141 


0.062 


0.013 


-0.056 


-0.083 


-0.116 


-0.160 


-0.196 


-0.218 


-0.208 


-0.171 


-0.127 


-0.091 


-0.069 


-0.058 


-0.050 


Cp 


-0.009 


-0.019 


-0.024 


-0.033 


-0.055 


-0.118 


-0.298 


-0.393 


-0.400 


-0.399 


-0.390 


-0.392 


-0.342 


-0.316 


-0.270 


-0.198 


-0.114 


-0.056 


-0.025 


-0.019 


-0.019 


-0.019 


-0.014 


B1 Antisinging Foil Experimental Data at AOA=5.2 deg 
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Chord ={455 mm 


X/C Y/C 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.1241 


-0.124 


-0.124 


-0.124 


-0.124 


- 0.120 


-0.115 


- 0.111 


-0.106 


- 0.102 


-0.097 


U/Uinf| 


1.001 


1.005 


0.995| 


0.989 


0.978 


0.987| 


1.003 


1.006 


1.006 


1.006| 


1.005 


Uinf= 5 


V/Uinf 


-0.049 


-0.046 


-0.046 


-0.042 


-0.029 


-0.008 


- 0.011 


-0.015 


-0.025 


-0.034 


-0.042 


.256 m/sec 


Cp 


- 0.002 


-0.006 


0.004 


0.010 


0.021 


0.012 


-0.003 


-0.007 


-0.007 


-0.007 


-0.005 


-0.093 

0.997 

-0.044 

0.002 

-0.088 

0.999 

-0.036 

0.000 

-0.084 

0.984 

-0.067 

0.014 

-0.080 

0.965 

-0.063 

0.033 

-0.075 

0.942 

-0.067 

0.054 

-0.071 

0.889 

-0.042 

0.104 

-0.071 

0.911 

0.044 

0.084 

-0.071 

0.957 

0.048 

0.041 

-0.071 

0.978 

0.042 

0.021 

-0.071 

0.986 

0.032 

0.014 

-0.071 

0.991 

0.023 

0.008 

-0.027 

0.997 

0.029 

0.003 

0.017 

| 0.999 

0.034 

0.001 


AOA=|2 


yj c 


-0.780 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


15 degrees 


Y/C U/Uinf 


0.061 


0.105 


0.105 


0.105 


0.105 


0.105 


0.105 


0.101 


0.096 


0.092 


0.087 


0.083 


0.078 


0.074 


0.069 


0.065 


0.060 


0.056 


0.052 


0.052 


0.052 


0.052 


0.052 


1.005 


1.008 


1.010 


1.010 


1.012 


1.029 


1.100 


1.183 


1.221 


1.246 


1.258 


1.265 


1.252 


1.240 


1.216 


1.161 


1.088 


1.027 


1.005 


1.008 


1.012 


1.016 


1.008 


V/Uinf 


0.032 


0.032 


0.046 


0.065 


0.093 


0.141 


0.185 


0.156 


0.108 


0.074 


0.008 


-0.017 


-0.053 


- 0.101 


-0.146 


-0.183 


-0.192 


-0.156 


-0.105 


-0.070 


-0.051 


-0.046 


-0.036 


Cp 


-0.005 


-0.009 


- 0.011 


- 0.012 


-0.017 


-0.040 


- 0.122 


- 0.212 


-0.252 


-0.279 


-0.291 


-0.301 


-0.285 


-0.274 


-0.250 


-0.190 


- 0.111 


-0.040 


- 0.010 


- 0.011 


-0.014 


-0.017 


-0.009 


B1 Antisinging Foil Experimental Data at AOA=2.15 deg 
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Chord = 


455 mm 




0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 










-0.082 


-0.0791 


-0.075 


-0.075 


-0.075 


-0.075 


-0.075 


-0.075 


-0.031 


0.0131 










0.999 


0.986 


0.937 


0.929 


0.962 


0.978 


0.989 


0.993 


0.997 


0.999 


Uinf= 5 
































-0.043 


-0.051 


-0.055 


0.018 


0.027 


0.023 


0.014 


0.010 


0.018 


0.019 


.134 m/sec 


1.35 degrees 


Y/C U/Uin 



B1 Antisinging Foil Experimental Data at AOA=1.35 deg 
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Chord =! 


XJ C 


455 mm 


Y/C 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.095 


-0.095 


-0.095 


-0.095 


-0.095 


-0.095 


-0.094 


-0.0931 


-0.092 


-0.091 


-0.090 


-0.089 


-0.088 


-0.087 


- 0.0861 


-0.085 


-0.084 


-0.084 


-0.084 


-0.084 


-0.084 


-0.084 


-0.040 


0.004 


U/Uinfl 


1.008 


1.004| 


1.006 


0.996 


0.986 


0.996 


1.019 


1.027 


1.037 


1.041 


1.049 


1.045 


1.054 


1.049 


1.049 


1.049 


1.016 


0.965 


0.977 


0.988 


0.992 


0.994 


0.998 


1.002 


Uinf=j5. 


V/Uinf 


-0.033 


-0.035 


-0.0291 


-0.027 


0.002 


0.031 


0.033 


0.025| 


0.021 


0.016| 


0.006 


0.004 


0.006 


-0.025 


-0.025 


-0.045 


-0.087 


-0.037 


-0.008 


0.000 


0.002 


0.002 


0.006 


0.006 


15 m/sec 


Cp 


-0.008 


-0.005 


-0.006 


0.004 


0.013 


0.003 


- 0.020 


-0.028 


-0.038 


-0.042 


-0.050 


-0.046 


-0.056 


-0.050 


-0.050 


-0.051 


-0.019 


0.034 


0.023 


0.012 


0.008 


0.006 


0.002 


- 0.002 


AOA=|-0. 


X/C 


1 degrees_ 


Y/C U/Uinf 


-0.780 


-0.780| 


-0.692 


-0.604 


-0.516 


-0.429| 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


0.048 


0.091 


0.091 


0.091 


0.091 


0.091 


0.091 


0.091 


0.090| 


0.089 


0.088 


0.087 


0.086 


0.085 


0.084 


0.083 


0.082 


0.081 


0.080 


0.080 


0.080 


0.080 


0.080 


1.002 


1.004 


1.002 


0.996 


0.988 


0.977 


1.004 


1.074 


1.188 


1.198 


1.202 


1.204 


1.190 


1.146 


1.082 


1.025 


1.004 


1.012 


1.014 


1.016 


1.012 


V/Uinf 


0.006 


0.010 


0.014 


0.021 


0.039 


0.072 


0.136 


0.153 


0.132 


0.107 


0.054 


0.033| 


-0.004 


-0.047 


-0.095 


-0.138 


-0.151 


-0.124 


-0.074 


-0.045 


-0.029 


-0.025 


- 0.021 


Cp 


- 0.002 


-0.004 


- 0.002 


0.004 


0.011 


0.020 


-0.013 


-0.088 


-0.141 


-0.182 


-0.208 


-0.218 


- 0.222 


-0.226 


-0.213 


-0.166 


-0.096 


-0.033 


-0.007 


-0.013 


-0.014 


-0.016 


- 0.012 


B1 Antisinging Foil Experimental Data at AOA=-O.l deg 










































































































































































































































































yjc\ 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


0.011 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


Y/C 


-0.079 


-0.079 


-0.079 


-0.079 


-0.079 


-0.080 


-0.081 


-0.082 


-0.083 


-0.084 


-0.086 


-0.087 


-0.088 


-0.089 


-0.090 


-0.091 


-0.092 


-0.092 


-0.092 


-0.092 


-0.092 


-0.092 


-0.048 


-0.004 


U/Uinfl V/Uinfl 


1.009 


1.013 


1.007 


1.000 


0.990 


1.000 


1.024 


1.039 


1.051 


1.060 


1.075 


1.073 


1.086 


1.096 


1.098 


1.113 


1.088 


0.994 


0.992 


0.996 


0.998 


1.000 


1.002 


1.002 


-0.023 


-0.025 


-0.023 


-0.015 


0.006 


0.041 


0.049 


0.043 


0.041 


0.038 


0.023 


0.021 


0.025 


- 0.011 


-0.013 


-0.045 


-0.126 


-0.094 


-0.051 


-0.030 


- 0.021 


-0.017 


-0.013 


-0.011 


CpI 


-0.009 


-0.013 


-0.007 


0.000 


0.010 


0.000 


-0.026 


-0.041 


-0.053 


-0.0621 


-0.078 


-0.076 


-0.090 


- 0.100 


-0.103 


- 0.120 


- 0.100 


0.002 


0.007 


0.004 


0.002 


0.000 


- 0.002 


- 0.002 


WC 


-0.780 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


Y/C 


0.040 


0.084 


0.084 


0.084 


0.084 


0.084 


0.084 


0.085 


0.086 


0.087 


0.088 


0.089 


0.090 


0.091 


0.092 


0.094 


0.095 


0.096 


0.097 


0.097 


0.097 


0.097 


0.097 


U/Uinfl V/Uinfl 


0.998 


1.000 


0.996 


0.988 


0.975 


0.952 


0.954 


1.020 


1.079 


1.1241 


1.152 


1.169 


1.183 


1.186 


1.175 


1.141 


1.0791 


1.022 


1.003 


1.007 


1.013 


1.013 


1.009 


- 0.011 


-0.009 


-0.008 


-0.006 


0.004 


0.034 


0.111 


0.149 


0.143 


0.124 


0.074 


0.057 


0.017 


-0.025 


-0.075 


- 0.121 


-0.141 


-0.117 


-0.070 


-0.043 


-0.030 


- 0.021 


- 0.021 


B1 Antisinging Foil Experimental Data at AOA=-l:5 deg 


Cp 


0.002 


0.000 


0.004 


0.012 


0.025 


0.046 


0.038 


-0.032 


-0.092 


-0.140 


-0.167 


-0.185 


-0.199 


-0.204 


-0.193 


-0.158 


-0.092 


-0.029 


-0.006 


-0.008 


-0.013 


-0.013 


-0.009 
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Appendix G 


B1 Cupped Trailing Edge Foil 
Experimental Data 




B1 Cupped Foil Geometry at AOA=0.0 
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Chord =! 


yjc 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.780 


455 mm 


Y/C 


-0.116 


-0.116 


- 0.1161 


-0.116 


-0.116 


- 0.112 


-0.109 


-0.106 


- 0 . 102 | 


-0.099 


-0.095 


-0.092 


-0.0851 


-0.081 


-0.078 


-0.075| 


-0.075 


-0.075 


-0.075 


-0.075 


-0.0751 


-0.031 


0.013 


0.057 


A0A=I2 


X/C 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


25 degrees 


Y/C U/Uinf 


U/Uinf 


0.997 


1.003 


0.999 


0.985 


0.975 


0.963| 


0.952 


0.961 


0.969 


0.975 


0.977 


0.973 


0.963 


0.948 


0.928| 


0.877 


0.908 


0.954 


0.979 


0.989 


0.997 


0.9991 


1.005 


1.011 


V/Uinf| 


-0.057 


-0.061 


-0.0631 


-0.061 


-0.049 


-0.051 


-0.026 


-0.016 


-0.008 


- 0.012 


- 0.020 


- 0.0221 


-0.041 


-0.039 


-0.041 


-0.016 


0.067 


0.069 


0.053 


0.043 


0.033 


0.037 


0.041 


0.041 


Cp 


0.002 


-0.005 


- 0.001 


0.013 


0.023 


0.035 


0.047 


0.038 


0.030 


0.024 


0.022 


0.026 


0.035 


0.050 


0.069 


0.115 


0.085 


0.043 


0.019 


0.010 


0.003 


0.000 


-0.006 


- 0.012 



0.094| 


0.091 


0.087| 


0.084 


0.081 


0.077 


0.074 


0.070 


0.067 


0.063 


0.060| 


0.060 


0.890 

0.060 

0.978 

0.060 

1.066 

0.060 


1.017 


1.018 


1.020 


1.022 


1.040 


1.115 


1.203 


1.247 


1.272 


1.272 


1.280 


1.270 


1.266 


1.247 


1.199 


1.125J 


1.052 


1.015 


1.005 


1.007 


1.007| 


1.001 


V/Uinf| 


0.041 


0.055 


0.075 


0.106 


0.157 


0.203 


0.173 


0.126 


0.090 


0.028 


0.002 


-0.043 


-0.090 


-0.144 


-0.185 


-0.208 


-0.191 


-0.149 


-0.106 


-0.083 


-0.067 


-0.059 


Cp 


-0.018 


- 0.020 


-0.023 


-0.028 


-0.053 


-0.142 


-0.239 


-0.285 


-0.313 


-0.310 


-0.319 


-0.308 


-0.306 


-0.287 


-0.236 


-0.154 


-0.071 


-0.026 


- 0.010 


- 0.010 


-0.009 


-0.003 


B1 Cupped Foil Experimental Data at AOA=2.25 deg 

































































































































































































































































B1 Cupped Foil Experimental Data at AOA=1.3 deg 































































































































































































































































Chord = 


Y1C 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.780 


455 mm 


Y/C| 


-0.095 


-0.095 


-0.095 


-0.095 


-0.095 


-0.094 


-0.093 


-0.092 


-0.091 


-0.091 


-0.090 


-0.089 


-0.087 


-0.086 


-0.085 


-0.085 


-0.085 


-0.085 


-0.085 


-0.085 


-0.085 


-0.041 


0.003 


0.047 


AOA=i 


X/C 


0.5 degrees 


Y/C 


U/Uinf 


0.991 


0.997 


0.989| 


0.981 


0.968 


0.956| 


0.951 


0.974 


0.983 


0.995 


1.004 


1.002 


1.006 


1.004 


1.002 


0.958 


0.935 


0.962 


0.981 


0.989 


0.993 


0.997 


1.000 


1.004 


V/Uinf 


-0.055 


-0.063 


-0.065 


-0.061 


-0.057 


-0.033 


-0.019 


0.002 


0.010 


0.000 


-0.019 


- 0.010 


-0.033 


-0.036 


-0.057 


-0.075 


- 0.002 


-0.004 


0.011 


0.015 


0.015 


0.011 


0.021 


0.021 


Cp 


0.008 


0.001 


0.009 


0.017 


0.030 


0.042 


0.048 


0.026 


0.017 


0.005 


-0.004 


- 0.002 


-0.007 


-0.005 


-0.004 


0.038 


0.063 


0.037 


0.019 


0.011 


0.007 


0.003 


- 0.001 


-0.004 


U/Uinf V/Uinf 


-0.780 

0.091 

1.006 

0.023 

-0.692 

0.091 

1.006 

0.031 

-0.604 

0.091 

1.002 

0.044 

-0.516 

0.091 

1.000 

0.065 

-0.429 

0.091 

1.000 

0.109 

-0.341 

0.091 

1.044 

0.170 

-0.253 

0.090 

1.123 

0.172 

-0.165 

0.090 

1.176 

0.142 

-0.077 

0.089 

1.211 

6.115 

0.011 

0.088 

1.236 

0.055 

0.099 

0.087 

1.251 

0.034 

0.187 

0.086 

1.251 

-0.006 

0.275 

0.085 

1.255 

-0.054 

0.363 

0.084 

1.243 

-0.107 

0.451 

0.083 

1.205 

-0.155 

0.538 

0.083 

1.144 

-0.182 

0.626 

0.082 

1.071 

-0.172 

0.714 

0.081 

1.029 

-0.132 

0.802 

0.081 

1.023 

-0.094 

0.890 

0.081 

1.023 

-0.069 

0.978 

0.081 

1.023 

-0.057 

1.066 

0.081 

1.018 

-0.050 


Cp 


-0.006 


-0.007 


-0.003 


- 0.002 


-0.006 


-0.060 


-0.145 


- 0.202 


-0.240 


-0.265 


-0.283 


-0.282 


-0.289 


-0.279 


-0.238 


-0.171 


-0.089 


-0.038 


-0.028 


-0.026 


-0.025 


-0.019 


B1 Cupped Foil Experimental Data at AOA=0.5 deg 
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Chord = 


x/c 


1.066 


0.978 


0.890 


0.802 


0.714 


0.626 


0.538 


0.451 


0.363 


0.275 


0.187 


0.099 


-0.077 


-0.165 


-0.253 


-0.341 


-0.429 


-0.516 


-0.604 


-0.692 


-0.780 


-0.780 


-0.780 


-0.780 


455 mm 


Y/C 


-0.0541 


-0.054 


-0.054 


-0.054 


-0.054 


-0.058 


-0.063 


-0.067 


-0.071 


-0.075 


-0.079 


-0.083 


-0.092 


-0.096 


- 0.1001 


-0.104 


-0.104 


-0.104 


-0.104 


-0.1041 


-0.104 


-0.060 


-0.016 


0.028 


U/Uin 


1.0201 


1.010 


1.004 


0.991 


0.981 


0.9721 


0.983 


1.008 


1.0291 


1.047 


1.070 


1.077 


1.110 


1.127 


1.162 


1.141 


1.0161 


1.000 


1.000 


1.002 


1.002 


1.000 


1.000 


0.998 


Uinf= 


V/Uinf 


-0.037 


-0.033 


-0.037 


-0.029 


-0.013 


0.012 


0.050 


0.065 


0.067 


0.069 


0.044 


0.058 


0.031 


0.019 


-0.019 


-0.140 


- 0.121 


-0.069 


-0.046 


-0.029 


-0.023 


-0.019 


-0.015 


-0.017 


5.198 m/sec 


Cp 


- 0.020 


- 0.011 


-0.005 


0.009 


0.019 


0.028 


0.016 


- 0.010 


-0.032 


-0.050 


-0.073 


-0.082 


-0.117 


-0.136 


-0.175 


-0.161 


-0.023 


-0.003 


- 0.001 


-0.003 


-0.003 


- 0.001 


- 0.001 


0.001 


AOA= 


yjc 


-0.780 


-0.692 


-0.604 


-0.516 


-0.429 


-0.341 


-0.253 


-0.165 


-0.077 


0.011 


0.099 


0.187 


0.275 


0.363 


0.451 


0.538 


0.626 


0.714 


0.802 


0.890 


0.978 


1.066 


-2.75 degrees 


Y/C U/Uinfl 


0.072 


0.072 


0.072 


0.072 


0.072 


0.072 


0.076 


0.080 


0.084 


0.088 


0.093 


0.097 


0.101 


0.105 


0.109 


0.113 


0.117 


0.122 


0.122 


0.122 


0.122 


0.122 


0.998 


0.9931 


0.985 


0.964 


0.925 


0.906 


0.979 


1.043 


1.093 


V/Uirtfl 


- 0.012 


-0.015 


-0.0171 


-0.013 


0.010 


0.096 


0.160 


0.1651 


0.154 


1.131 

0.112 

1.158 

0.096 

1.179 

0.058 

1.197 

0.013 

1.200 

-0.040 

1.175 

-0.096 

1.127 

-0.137 

1.060 

-0.135 

1.025 

-0.094 

1.018 

-0.062 

1.020 

-0.042 

1.020 

-0.035 

1.014 

-0.031 


Cp 


0.001 


0.007 


0.015 


0.035 


0.072 


0.085 


0.008 


-0.057 


-0.109 


-0.146 


-0.175 


-0.197 


-0.216 


- 0.221 


-0.195 


-0.145 


-0.071 


-0.030 


- 0.020 


- 0.021 


- 0.020 


-0.014 



B1 Cupped Foil Experimental Data at AOA=-2.75 deg 
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Appendix H 

Decaying Vortex Tracking Program 
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